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Vorwort 



Die vorliegende Arbeit ist ein Forschungsbericht über das 
vom Ministerium für Wissenschaft und Forschung des Landes 
Nordrhein-Westfalen in den Jahren 1975 - 77 unterstützte 
Projekt "Konstruktion und Einsatz von Digitalfiltern zur 
Analyse und Prognose ökonomischer Zeitreihen". 

Der Bericht ist bev,nißt knapp gehalten, da zum einen eine 
vollständige Darstellung aller erarbeiteten Resultate den 
vom Ministerium gesteckten Umfang bei weitem überschrit- 
ten hätte und zum anderen eine solche ausführliche Dokumen- 
tation an dieser Stelle nicht erforderlich ist, da über 
die wichtigsten Aspekte des Projektes separate Publika- 
tionen erfolgen. Außerdem waren manche Verfahren beim 
offiziellen Abschluß des Projektes noch nicht anv/endungs- 
reif. 

Bisher sind vom Proj ektleiter folgende projektbezogene 
Arbeiten erschienen bzw. befinden sich im Druck: 

Zur Diskussion der sog. "langen Konjunkturwellen" - eini- 
ge statistische Bemerkungen, in: Wirtschaftsdienst (De- 
zember 1 976) . 

Grundprobleme der Saisonbereinigung ökonomischer Zeitrei- 
hen, in: Jahrbücher für Nationalökonomie und Statistik, 

Bd. 1 92, Heft 3-4 (1977) . 

Über eine Klasse von einfachen FIR-Tief paß-Selektions- 
Filtern, in: Allgemeines statistisches Archiv, Heft 3, 

1 978. 

Zur Verknüpfung von "Zeitbereichsanalyse" und "Frequenz- 
bereichsanalyse" - ein Beitrag zur Zeitreihenanalyse, in; 
Mathematische Operationsforschung und Statistik, Ser. 
STATISTICS 9 (1978) Nr. 3. 




IV 



In Vorbereitung befinden sich eine größere Arbeit von H. 
Schulte zum Problem der "langen Konjunkturwellen" ("Metho- 
dische Untersuchungen zum Problem langer Wellen"), die bis 
zum Sommer dieses Jahres abgeschlossen sein wird, sowie 
zwei Arbeiten des Projektleiters , welche die im vorliegen- 
den Bericht skizzierten Verfahren zur Saisonbereinigung 
ökonomischer Zeitreihen ausführlich darstellt. 

Mein besonderer Dank gilt Herrn Dipl. -Math. Heinrich Schulte, 
ohne dessen kritische Mitarbeit und engagiertem Einsatz das 
Forschungsprojekt nicht durchzuführen gewesen wäre. 

Zu Dank verpflichtet bin ich auch Frau Ruth Eckert für die 
sorgfältige Erstellung eines reproduktionsfähigen Manuskrip- 
tes sowie Herrn Kollegen Prof. Dr . Arno Jaeger für seine 
großzügige Überlassung von "Schreibkapazität". 

Bochum, im April 1978 Winfried Stier 
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1 . Nicht-rekursive Filter 



1.1. Grundbegriffe der Filtertheorie 



Unter einem "Filter" versteht man einen Algorithmus zur Trans- 
formation einer Zeitreihe ("Signal") : 



y(t) = F[x(t) ] , t=1 ,2, . . .T 



Dabei ist x(t) die vorliegende Zeitreihe, y(t) die gefilterte 
Zeitreihe und T die Länge von x (t) . Prinzipiell ist es dabei 
gleichgültig, ob dieses Signal als deterministisch oder als Re- 
alisation eines stochastischen Prozesses angesehen v;ird. Den 
theoretischen Hintergrund der hier darzustellenden Filter lie- 
fert die Theorie der diskreten, linearen, zeitinvarianten 
Systeme. (Vgl. dazu etwa FREEMAN , 9) Ein in der Systemtheorie 
v/ichtiges Testsignal ist die sog. "diskrete Impulsfunktion": 



6 (t) 



. 0 für t 40 
^ 1 für t=0 



Damit läßt sich jedes Signal x(t) in der Form darstellen: 
X (t) = Ex (k) 6 ( t-k) 

k=~oo 

Für das gefilterte Signal ergibt sich dann: 
y(t) = F[ E x(k)ö(t-k)] 

k=-oo 



= ■ Ex(k)F[6 (t-k) ] 

k=-oo 



= E X (k) h (t-k) 

k=-oo 



= Eh (k) X (t-k) 
k=-cx) 



Dabei ist 

h(t) :=F[6 (t) ] 

die sog. "Impuls-Antv 7 ortf unktion" des Filters F. Durch h(t) ist 
also der Filter F vollständig charakterisiert. Ist h(t) von end 
lieber Dauer , dann spricht man von einem "Finite-Impulse-Respon 
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se”-Filter (FIR-Filter ) , bei nicht endlicher Dauer dagegen von 
einem " Inf inite-Impulse-Response" -Filter (IIR-Filter) oder auch 
kurz von einem "rekursiven Filter" (Vgl. RABINER et al. / 31) . 

Für Filter, die auf der Theorie der oben charakterisierten Sy- 
steme basieren, kann die sog. "Filtergleichung" (das ist die 
Relation zwischen Filter-Input und Filter-Output) in Form einer 
linearen Differenzengleichung mit konstanten Koeffizienten an- 
gegeben werden : 

M N 

^ ^ ^ x(t-r) , N=t=M oder N=M 

k=0 r=0 ^ 



mit x(t);= Filter-Input 
y(t);= Filter-Output 



zum Zeitpunkt t 
zum Zeitpunkt t 



Mit Hilfe der z-Transf ormation : 

X(z):=Z[x(t)]= ? x(t)z ^ , z e C 

t=-oo 

kann für obige Differenzengleichung geschrieben werden: 



^ -k 

E a, z ^Y(z) 

k=0 



E b z ^X(z) 
r 

r=0 



Die Funktion 



T(z) 



Y(z) 

X(z) 



E b^z 
r=0 
M 

E a, z 
k=0 



1 



ist für Werte auf dem Einheitskreis, also für z=e^^, die sog. 
"Frequenz-Antwortfunktion" des obigen Filters. T(e^^) ist im 
allgemeinen eine komplexe Funktion mit der Periode 2tt und ist 
deshalb in der Form 



T(e^^ = 






|A|ä^ 
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darstellbar, wobei |T(e^^) | der Betrag und cp(X) die Phasenfunk- 
tion des Filters ist. An der Funktion T(e^^) (kurz: T(X)) lassen 
sich relativ bequem die Ausblendeigenschaften eines konkreten 
Filters studieren. Umgekehrt können aber auch Filter so konstru- 
iert werden, daß sie eine vorgegebene Frequenz -Antwort funktion 
bzw. Transferfunktion (darunter wird das Betragsquadrat von T(X) 
verstanden ) "optimal" approximieren. Zur Lösung der anste- 
henden Approximationsprobleme können verschiedene V'Jege einge- 
schlagen werden, die sich u. a. auch dadurch unterscheiden, ob 
eine vorgegebene Transferfunktion mit Hilfe eines IIR-Filters, 
also rekursiv, oder mit Hilfe eines FIR-Filters, also nicht- 
rekursiv, realisiert werden soll. Obv/ohl in vorliegender Ar- 
beit rekursive Filter im Mittelpunkt stehen, soll auf FIR-Fil- 
ter zu Beginn der Arbeit eingegangen werden, wobei zwei Kon- 
struktionsprinzipien erläutert v/erden und die Frage untersucht 
wird, ob derartige Filter praktisch und theoretisch einen Ge- 
winn für die Analyse von ökonomischen Zeitreihen dar stellen. 

1.2. Konstruktion von FIR-Filtern durch 
"windowing " 



Als praktisch wichtig erv/eist sich die Konstruktion von sog. 
"Selektionsfiltern". Das sind Filter, deren Frequenz -Antwort- 
funktion in einem gewissen Frequenzintervall (oder in mehre- 
ren Intervallen) den V7ert Eins einnehmen und sonst den Wert 
2 ) 

Null . Auf solche Filter führen häufig substanzwissenschaft- 
lich motivierte Problemstellungen, wie sie aus der Zeitreihen- 
analyse bekannt sind. Im Rahmen solcher Analysen tritt z. B. 
das Problem auf, gewisse "Komponenten" einer konkreten Zeit- 
reihe zu schätzen. Dabei erfolgt die jeweilige inhaltliche 
Deutung dessen, was unter einer "Komponente" verstanden wer- 
den soll bzv7. v/ie mehrere Komponenten gegeneinander adäquat 
abzugrenzen sind im Rahmen substanzwissenschaftlicher Erör- 
terungen. So ist* etwa die Schätzung der sog. "glatten Kompo- 
nente" eine bekannte Problemstellung bei der Analyse von 

1) Die beiden Begriffe werden in der Literatur nicht immer 
streng unterschieden, da es oftmals gleichgültig ist, wel- 
cher Begriff zugrunde gelegt wird. Auch in dieser Arbeit 
wird gelegentlich so verfahren. 

2) Oder allgemeiner: Filter, deren Frequenz -Antv/ort bzv/. 
Transf erf unkt ion intervallweise konstant sind. 
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wirtschaftlichen Zeitreihen, wobei von der Vorstellung ausge- 
gangen wird, daß diese Komponente die in der ökonomischen Rea- 
lität sich vollziehenden Wachstums- und Konjunkturprozesse 
reflektiert. Leider reichen jedoch in der Regel derartige sub- 
stanzwissenschaftliche Überlegungen nicht aus, um zu prakti- 
kablen Schätzverfahren zu gelangen, da sie meistens zu "un- 
scharf" sind. Daher ist der Statistiker gezwungen, mit forma- 
len Kriterien zu arbeiten, die es ihm gestatten, operable 
Komponenten-Def initionen einzuführen, die zu realisierbaren 
Schätzverfahren führen, wobei natürlich zu fordern ist, daß 
diese Definitionen auch ökonomisch sinnvoll sind und sich 
möglichst mit dem decken, was dem Wirtschaftstheoretiker vor- 
schwebt . 

Da es sich auch für wirtschaftliche Zeitreihen als zweckmäßig 
erwiesen hat, sie als Realisationen von stochastischen Prozes- 
sen zu betrachten und die Analysen nicht nur im Zeitbereich, 
sondern auch im Frequenzbereich durchzuführen, erscheint für 
die in vorliegender Abhandlung interessierende Problemstellung 
"Schätzen der glatten Komponente einer ökonomischen Zeitreihe" 
als formales Definitionskriterium der relevanten Komponente 
eine adäquate Zerlegung des Frequenzintervalls [0,tt] in zwei 
diskunkte Intervalle = [ O, A* ] , 1^= (X* tt ] sinnvoll. Dabei sind 
alle Komponenten, die dem Intervall (A*,tt] zuzuordnen sind, 
als "nicht-glatt" zu bezeichnen. Allerdings ist eine solche 
Zerlegung nicht eindeutig, da a priori nicht feststeht, wie 
A* zu wählen ist. Aus ökonomischen Gründen muß A* größer als 

Null sein, da die glatte Komponente auch konjunkturell beding- 

2 TT 

te Schwingungen umfassen soll und kleiner als jj' saiso- 
nale Schwankungen in einer Zeitreihe ex definitione nicht zu 
den langfristigen, durch Wachstums- und Konjunkturprozesse 
hervorgerufenen Bewegungen einer wirtschaftlichen Zeitreihe zu 
zählen sind. Um zu einer eindeutigen Definition der glatten 

Komponente im Frequenzbereich zu gelangen, muß also jeweils 
2tt 

ein A*e(0,y2) festgelegt werden. Anstatt die glatte Komponen- 
te im Frequenzbereich zu definieren, kann jedoch auch eine 
Definition im Zeitbereich vorgenommen werden. Dabei wird üb- 
licherweise angenommen, daß diese Komponente innerhalb eines 
Bereiches, dessen Länge vorzugeben ist, durch Funktionen aus 




5 



einem linearen Raum approximiert werden kann, wobei vorzugs- 
weise Polynome Verwendung finden, deren Grad ebenfalls vor- 
zugeben ist. Damit ist dann eindeutig definiert, was unter 
einer glatten Komponente zu verstehen ist. Bei diesem Ansatz 
werden somit ganz bestimmte Tra j ektorieneigenschaf ten für die 
glatte Komponente einer wirtschaftlichen Zeitreihe unterstellt, 
die letzten Endes auf die Hypothese hinauslaufen, man könne 
diese Trajektorie durch einen geschlossenen mathematischen 
Ausdruck darstellen. Dies ist sicher eine sehr restriktive 
Hypothese. Allerdings hat sie den Vorzug, daß die so defi- 
nierte glatte Komponente durch eine lineare erwartungstreue 
Schätzfunktion mit minimaler Varianz ermittelt werden kann, 
wobei diese Schätzfunktion relativ einfach in Form von glei- 
tenden Durchschnitten darstellbar ist. (Für Einzelheiten 
siehe WEICHSELBERGER, 40) 

Ist man nun aber der Auffassung, daß eine derartige Hypo- 
these für eine vorliegende Zeitreihe unzweckmäßig oder zu 
eng ist, dann muß ein anderer Ansatz gewählt werden. Im fol- 
genden soll ein entsprechendes Schätzverfahren diskutiert 
werden, dessen theoretische Basis im Frequenzbereich formu- 
liert ist. An die Stelle von "Erwartungstreue" und "minima- 
ler Varianz" treten dabei andere Optimalitätskriter ien , die 
aus der Theorie des "f ilter-designs" stammen. 

Das eben spezifizierte Problem der Ermittlung der glatten 
Komponente einer ökonomischen Zeitreihe erfordert einen Se- 
lektionsfilter in Form eines "Tiefpasses". Für einen "idea- 
len" Tiefpaß gilt 

1 für 

T*(A) = { ^ . 

O für A^< 1 A I < TT ^ 

Dabei ist A^ die sog. "Abschneidefrequenz" des Filters 
( " cut-of f " -f requency) . Sie ist für die vorliegende Problem- 
stellung identisch mit der oben erwähnten Frequenz A’*'. 



1) Für diesen Filter sind Frequenz-Antwortfunktion und 
Transferfunktion identisch. 




6 



Die Konstruktion eines Tiefpaßfilters ist nun auf verschiede- 
ne Arten möglich. Gilt für die Koeffizienten im Nenner der 
Funktion H(z) 



^ r 1 für k=0 
k ^ 0 " k+0 

dann lautet die entsprechende Filtergleichung 

y(t) = b^x (t) +b^x (t-1 ) + . . . +b^x (t-N+1 ) 

d. h. der Filter-Output zum Zeitpunkt t hängt von den Werten 
des Filter-Inputs zum Zeitpunkt t und den vorangegangenen 
Zeitpunkten t-1 , . . . , t-N+1 ab. Ein solcher Filter heißt "nicht 
rekursiv" oder auch "Transversalfilter". Für einen "rekursi- 
ven" Filter dagegen lautet die Filtergleichung 



y(t) = a^ y (t-1) t . . .+aj^y (t-M) +b^x (t) . . .+bj^x (t-N+1 ) 



d. h. der Filter-Output zum Zeitpunkt t hängt auch vom Fil- 

1 ) 

ter-Output der Zeitpunkte t-1,...,t-M ab. (Vgl. RÄDER, 
32) Die hier zu diskutierende Klasse von Tiefpaßfilter sol- 
len als Transversalfilter implementiert werden. Das hat zur 
Folge, daß diese Filter FIR-Filter sein müssen, da ja in 
praxi N, d.h. also die Filterlänge, immer nur endlich sein 
kann . 



Für einen transversalen FIR-Filter der Länge N kann also 
geschrieben werden 

N-1 , 

T (A) = E b e , I A |<TT 

r=0 ^ 

und 

N-1 

y (t) = E b X (t-r) 
r=0 



1) Bei einem Transversalfilter hat H(z) offensichtlich nur 
Pole für z=0, bei einem rekursiven Filter existieren da- 
gegen auch Pole für z4o. 
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Andererseits gilt 

oo 

y(t) = E h(k)x(t-k) 
k=0 

wenn der Einfachheit halber von einem "kausalen" Filter aus- 
gegangen wird, für den gilt 

h(k) =0 für k<0 ^ ' . 

Dabei ist h(k) die Impuls-Antwortfunktion des zu konstruie- 
renden Tiefpaßfilters. Aus den letzten beiden Ausdrücken er- 
geben sich folgende Beziehungen zwischen den b^ und den h(k): 

b^ = h(0) , b^ = h(1) , . . . = h(N-1) 

b^ = h(N) = O , b^^^ = h(N+1) = 0,... 

Somit kann für die Filtergleichung geschrieben werden 

N-1 

y(t) = E h(k)x(t-k) . 
k=0 

Die N unbekannten Koeffizienten der Filtergleichung sind al- 
so identisch mit den N ersten Gliedern der Impuls-Antwort- 
funktion des Filters. Diese läßt sich relativ einfach mit Hil- 
fe der Theorie der Fourier-Reihen bestimmen: 

Da T(A) eine periodische Funktion in A ist mit der Periode 
2tt , kann sie formal in eine unendliche Fourier-Reihe ent- 
wickelt werden: 

oo 

T(A) = E a(k)e"^^^ , |X|<tt 

k=-oo 



1) Diese Annahme wird hier getroffen, da die Filtergleichung 
in nicht-symmetrischer Form vorliegt. Die nachfolgenden 
Ableitungen sind jedoch davon unabhängig. In den Simula- 
tionsstudien wird die Kausalitätsbedingung wieder aufge- 
geben . 

2) Zur Theorie der Fourier-Reihen vgl. etwa: KNOPP, 18. 
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Für die Fourier-Koeffizienten a (k) gilt: 

a(k) / T(A) k=0 , ±^ , ±2 , . . . 

-TT 

Da aber außerdem gilt: 



-irA 

T(A) = Z b e , |A|<tt 

r=0 



folgt mit den obigen Identitäten: 



h(0) = a(0) , h(1) = a(1) , . . . ,h(N-1) = a(N-1) 



Die ersten N Fourier-Koeffizienten von T(A) stellen somit 
die Impuls-Antwortfunktion des FIR-Filters dar. 

Für einen idealen Tiefpaß ergeben sich für die Filter-Koef- 
fizienten h(k): 

h(k) = f r ^ ^ 2, . . . 

‘^A 

sin A k 

_ ^ _ I 1 1 o 



sin 27Tf^k 

TTk 



, k=+1, +2,... 



mit f , : = 2tt A_ 

A A 

und h(0) = 2f^ 



Ein idealer Tiefpaßfilter ist somit eigentlich ein IRR-Fil- 
ter. Da aber hier ein Tiefpaßfilter als Transversalfilter 
implementiert werden soll, der nur als FIR-Filter realisier- 
bar ist, können nur N Koeffizienten h(k) berücksichtigt wer- 
den. Dies ist gleichbedeutend damit, daß die Transferfunktion 
dieses Filters nicht ideale Rechteckgestalt aufweist, da 
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ihre unendliche Fourier-Reihe nach N Gliedern abgebrochen wird. 
Aus der Theorie der Fourier-Reihen ist bekannt (Vgl. etwa 
BARY, 3 ) , daß ein solches Vorgehen äquivalent ist mit einer 
Approximation der idealen Transferfunktion im Sinne des fol- 
genden Optimalitätskriteriums : 



1 2 
^ / [T(A) -T* (A) ]^dA = Min! 

ZTT 

-TT 

Ein so konstruierter Tiefpaßfilter ist also "optimal" in 
dem Sinne, daß dieses im Frequenzbereich formulierte Opti- 
malitätskriterium des "minimalen mittleren quadratischen 
Approximationsfehlers" erfüllt ist. 

Die Ersetzung einer unendlichen Fourier-Reihe durch ihre Par- 
tialsumme der Länge N kann formal durch Einführung einer sog. 
"Gewichts"- oder "Fenster "-Funktion mit den Eigenschaften 



WR(k) 



1 für k = 0,1,2,...N-1 
0 für K >N-1 



dargestellt werden. Für die Filter-Koeffizienten läßt sich 
damit schreiben 



hj^(k) = h(k) Wj^(k) 



und der entsprechende Filter soll als "Rechteck"-Filter be- 
zeichnet werden. Die realisierte Transferfunktion läßt sich 
dann als Faltungsintegral darstellen 



T„(») 



^ ; T*(X)Wj^(A-0)de 

-TT 



wobei ( A- 0) 

K 



sin (A-e)N/2 
sin (A-e) /2 
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die Fouriertransf ormierte der Fensterfunktion Wj^(k) ist. 

1.2.1. Ein Simulationsbeispiel 

Die Wirkung eines solchen Tiefpaßfilters soll nun an einer si- 
mulierten Zeitreihe demonstriert werden. Diese Zeitreihe wur- 
de wie folgt erzeugt: 

X(t) = 100+0,001 1^-0, 1 1^+0 ,5t+200sin0,06TTt+e (t) 
mit e (t) ~ N(0;50) 

Die glatte Komponente dieser Reihe setzt sich also additiv 
aus einem Polynom dritter Ordnung und einer niederfrequen- 
ten Schwingung zusammen. Diesem "Signal", das hier als nicht- 
stochastisch angenommen wird, ist ein unkorrelierter Zufalls- 
prozeß ("weißes Rauschen") überlagert. 

2 ) 

In Abb. 1 sind sowohl die Originalreihe als auch die aus- 
gefilterte glatte Komponente wiedergegeben. Verwendet wur- 
de ein symmetrischer Tiefpaßfilter mit f^ = 0,05 und N = 23, 
so daß Phasenbetrachtungen außer acht bleiben können. Die 

Güte der Schätzergebnisse wurde mit Hilfe des Ungleichheits- 

3) 

koef f izienten von Theil gemessen . Bei diesem Filter fällt 
auf, daß er offensichtlich nicht in der Lage ist, den in 
der Reihe vorhandenen starken Trend zu erfassen. Der Trend 



1 ) Dies sind bekannte Zusammenhänge aus der Theorie der 
Fourier-Reihen. W_(A) ist dort als "Dirichlet-Kern " be- 

K. 

kannt . 



2 ) 

3 ) 



Alle Abbildungen dieses Berichtes befinden sich im Anhang. 



TH = 



[E(G -G 
t ^ ^ 

[EGj]'/^[Ea2]V2 
t ^ t ^ 



, 0<TH<1 



G^ = theoretische glatte Komponente von X(t) 
G^ = ausgefilterte glatte Komponente von X(t) 




wird systematisch unterschätzt und zwar zunehmend mit grö- 
ßer v/erdendem t. Der Grund für diesen "bias" kann leicht er- 
mittelt werden. Er ist in einer Eigenschaft der Transferfunk- 
tion des Rechteck-Filters an der Stelle A=0 zu suchen. 

Für die ideale Transferfunktion T*(A) gilt: 



T* (A) 



A=0 



E h(k) 



k=-oo 



1 



da T*(A) die Dirichletschen Bedingungen (zu diesen Bedingun- 
fen vgl. LANCZOS, 21 ) für die Konvergenz einer Fourier- 
Reihe erfüllt. Ein in der Reihe vorhandener Trend wird so- 
mit vollständig durch den idealen Tiefpaßfilter erfaßt. Beim 
Rechteckfilter dagegen ist dies im allgemeinen nicht mehr 
der Fall, da seiner Transferfunktion nur eine endliche Fourier- 
Reihe zugeordnet ist. Deshalb erweist es sich für die Analy- 
se von ökonomischen Zeitreihen als unumgänglich, die Filter- 
koeffizienten auf die Summe Eins zu normieren. Abb. 2 zeigt, 
daß damit eine systematische Unterschätzung des Trends ver- 
mieden wird. Auch mit Normierung weist der Rechteckfilter noch 
einige Eigenschaften auf, die u. U. für praktische Zwecke 
nachteilig sind: 

Die realisierte Transferfunktion stellt einen "inte- 

gralen Durchschnitt" aus der Fuktion Wj^(A) und der idealen 
Transferfunktion T*(A) dar. V1^{X) weist aber außer dem Haupt- 
maximum an der Stelle A = O noch mehrere Nebenmaxima bzw. 
Nebenminima (" side-lobes" ) auf . Das hat zur Folge, daß der 
Rechteckfilter möglicherweise Komponenten miterfaßt (wenn 
auch nur unvollständig) , die frequenzmäßig dem Sperrbereich 
^A^ I ^ 1 1"^ zuzuordnen sind ( "Filter-leakage " ) , da sich diese 
"side-lobes" auch in T (A) wiederfinden. Außerdem konvergiert 
die Transferfunktion des Rechteckfilters an den Stellen 
A = +A^ nicht gleichmäßig gegen die ideale Transferfunktion, 
da diese dort unstetig ist. Gemäß den Dirichletschen Bedin- 
gungen konvergiert sie dort gegen den Wert 0,5 und nicht ge- 
gen den Wert 1. Aufgrund dieser Unstetigkeiten läßt sich u. U. 




das aus der Theorie der Fourier-Reihen bekannte Gibbsche Phä- 
nomen beobachten (Vgl. dazu KNOPP, 18): T (X) verläuft oszil- 
latorisch um T*(X), wobei die Amplitude dieser Oszillation auch 
mit zunehmender Filterlänge N nicht abnimmt. In Abb. 3 ist dies 
für N = 99 und = 0,05 veranschaulicht . 

Tiefpaßfilter, welche diese Nachteile nicht oder nur sehr ab- 
geschwächt auf weisen, lassen sich nun auf zweierlei Arten er- 
zielen. Die Ursache für das Auftreten von ausgeprägten "side- 
lobes" bei T 2 (X) liegt nicht nur daran, daß diese primär bei 
W_ (X) auf treten, sondern auch daran, daß der Übergang vom 

X\ 

Durchlaß- zum Sperrbereich bei der idealen Transferfunktion 
T*(X) abrupt erfolgt.. Statt dessen könnte man aber auch einen 
"glatteren" Übergangsbereich definieren, etwa so, daß die 
ideale Transferfunktion nicht mehr Rechteck- sondern Trapez- 
gestalt aufweist. Dann ist die Transferfunktion nicht mehr 
unstetig an den Stellen X = + X^. Wichtige Filter dieses 
Typs sind die sog. "Ormsby "-Filter . (Vgl. ORMSBY, 26) Eine 
andere Möglichkeit besteht darin, statt eines Rechteckfensters 
andere Fensterfunktionen einzuführen. Wie aus der Theorie 
der Fourier-Reihen bekannt ist, hat das zur Folge, daß dadurch 
die Konvergenzeigenschaften der realisierten Transferfunk- 
tionen an den Unstetigkeitsstellen wesentlich verbessert wer- 
den. Diese Möglichkeit soll im folgenden betrachtet werden. 
Prinzipiell können hier alle Fenster Verwendung finden, deren 
Einsatz auch in der Spektralanalyse stationärer stochasti- 
scher Prozesse üblich ist. (Eine ausführliche Darstellung 
derartiger Fenster findet sich z. B. bei PARZEN, 28) 



.1) Wie man zeigen kann, treten diese Oszillationen nur dann 
auf, wenn bei gegebenem kleinen X^ die Filterlänge N ent- 
sprechend groß ist bzw. bei gegebener kleiner Filterlänge 
N die Abschneidefrequenz X^ entsprechend groß ist. 
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Für vorliegende Simulations Studie wurden folgende Fenster 
benutzt: 

1 . Bartlett-Fenster : 



WßCk) 



2 |k 
N-1 



|k| 



< 



N-1 

2 



2 ) 



2 . Hamming-Fenster ; 



w 

HAM 



(k) 



0,54+0,46 cos 



|k| 



< 



N-1 

2 



3 . Hanning-Fenste r i 



" 0,5+0, 5 cos , |k| < 

4 . Blackmann-Fenster : 

^BL ^ C,08 cos , |k| < 



5. Parzen-Fenstöi 



3) 



1 - 6 (irr) ^+6(^)3 



w (k) = { 



2 (1 -|M)3 



'N-1 
N-1 



0 <Jk| 






1) Dieses Fenster geht eigentlich auf Fejer zurück (Vgl. 
fej^;r, 7 ) 

2) Aus naheliegenden Gründen wird hier immer eine ungerade 
Filter länge gewählt. 

3) Dieses Fenster geht eigentlich auf Jackson und de la 
Vallee-Poussin zurück. (Vgl. ACHIESER, 1) 
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Ein weiteres Fenster, das bisher in der Spektralanalyse nicht 
verwendet wurde ist das 



6 



Kaiser-Fenster (Vgl. KAISER 

,(ß\l 



16) : 



Wj^(k) = 



1 - (— ^) ^ 



lo(3) 



< 



N-1 

2 



Dabei ist Iq(x) eine modifizierte Besselfunktion erster Art 
und nullter Ordnung. Dieses Fenster ist '‘optimal" in dem Sinn 
daß sein "Hauptfenster " ( "main-lobe" ) für eine vorgegebe- 

ne maximale "side-lobe" -Amplitude den größten Teil der Gesamt 
fläche des Fensters enthält. Dabei ist 2 < ß < 10 ein Parameter 
mit Hilfe dessen ein "trade-off" zwischen der "Breite" des 
Hauptfensters und der maximalen "side-lobe "-Amplitude vor- 
genommen werden kann . 

Die Abb. 4 zeigt beispielhaft die Originalreihe und mit Hil- 

2 ) 

fe des Kaiser-Filters geschätzte glatte Komponente 



Die Koeffizienten dieses Filters lauten: 



h( 


1) = 


-0.00001 = h(23) 


h( 7) = 0.03358 = h(17) 


h( 


2) = 


0.0 = h(22) 


h( 8) = 0.05632 = h(16) 


h( 


3) = 


0.00042 = h(21) 


h( 9) = 0.08270 = h(15) 


h( 


4) = 


0.00225 = h(20) 


h(10) = 0.10790 = h(14) 


h( 


5) = 


0.00714 = h(19) 


h(11) = 0.12616 = h(13) 


h( 


6) = 


0.01711 = h(18) 


h(12) = 0.13285 = h(12) 



1) Darunter ist der symmetrisch um k = 0 - liegende Teil des 
Fensters bis zu seiner ersten Nullstelle zu verstehen. 

2) Auf eine Wiedergabe der Grafiken der mit den anderen 
Fenstern erzielten Ergebnisse sei hier aus Platzgrün- 
den verzichtet, da sich diese optisch kaum von obiger 
Grafik unterscheiden. Wie beim Rechteckfilter ist auch bei 
den Fenster-Filtern eine Normierung der Filter-Koeffi- 
zienten auf die Summe Eins vorgenommen worden. 
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Die nachstehende Tabelle enthält die Theilschen Ungleichheits- 
koeffizienten für alle sieben "Fenster "-Filter : 

Filter Ungleichheitskoeffizienten 



von Theil 



Rechteck 


0,01311 


Bartlett 


0,00986 


Hamming 


0,00954 


Hanning 


0,00920 


Blackmann 


0,00826 


Parzen 


0,00787 


Kaiser (3=9) 


0,00817 



Offensichtlich werden durch die Einführung von nicht-recht- 
eckigen Fensterfunktionen beträchtlich bessere Schätzergeb- 
nisse gegenüber dem Rechteckfilter erzielt. Weitere hier nicht 
wiedergegebene Simulationsstudien zeigen, daß wie in der Spek- 
tralanalyse von einer "besten" Fensterfunktion im allgemeinen 
nicht gesprochen werden kann. 



1.2.2. Ein praktisches Beispiel 

Schließlich sei hier noch der Einsatz von FIR-Filtern an 
einem praktischen Beispiel, dem "Index der industriellen Net- 
toproduktion, gesamte Industrie" (Dezember 1962 - Januar 1975), 
demonstriert. Verwendet man z. B. einen Parzen-Filter 
(N = 23, f^ = 0,05) mit den Gewichten: 



h( 


1) 


= 0,0 


= h 


(23) 


h( 


7) 


= 0, 


,02924 = 


h 


(17) 


h( 


2) 


= 0,0 


= h 


(22) 


h( 


8) 


= 0, 


,05316 = 


h 


(16) 


h( 


3) 


= 0,00019 


= h 


(21) 


h( 


9) 


= 0, 


,08225 = 


h 


(15) 


h( 


4) 


= 0,00135 


= h 


(20) 


h(10) 


= 0, 


,11117 = 


h 


(14) 


h( 


5) 


= 0,00502 


= h 


(19) 


h(11) 


= 0, 


,13325 = 


h 


(13) 


h( 


6) 


= 0,01344 


= h 


(18) 


h(' 


12) 


= 0, 


,14186 = 


h 


(12) 



so fällt bei der geschätzten glatten Komponente dieser Reihe 
auf, daß diese einen, wenn auch nur schwach ausgeprägten zy- 
klischen Verlauf zeigt (s. Abb. 5) . Der Grund dafür ist darin 
zu suchen, daß die Transfer funktion dieses FIR-Filters bei der 
ersten saisonalen Frequenz f = 0,083 hinsichtlich ihrer Dämp- 
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fungseigenschaf ten unbefriedigend verläuft (Vgl. Abb. 6) . 
Dieser Filter ist deshalb für Monatsdaten mit ausgeprägter 
saisonaler Variation nicht verwendbar, dagegen ohne weiteres 
für Quartals-, Halbjahres- und Jahresdaten. 

Der Filter könnte nun dadurch verbessert werden, daß eine grö 
ßere Filterlänge gewählt wird. Da aber damit unvermeidlich 
weitere Reihenwerte verlorengehen, ist ein solches Verfahren 
wohl kaum zu empfehlen. Zweckmäßiger ist es, einen anderen 
Filter zu verwenden, etwa einen aus der Klasse der Kaiser- 
Filter. Da bei ihnen über einen weiteren Parameter verfügt 
werden kann, nämlich 3, zeichnen sich diese Filter durch eine 
gewisse Flexibilität aus. Wählt man z. B. 3 2 (aber 3 > 2) , 



dann 


ergeben sich 


die Gewichte 


(wieder 


mit 


N = 23 


' "a = ^ 


h( 


1) 


= 


-0.004 


= h(23) 


h( 7) 


= 


0.052 


= h(17) 


h( 


2) 


= 


0.0 


= h(22) 


h( 8) 


= 


0.065 


= h(16) 


h( 


3) 


= 


0.006 


= h(21) 


h( 9) 


= 


0.077 


= h(15) 


h( 


4) 


= 


0.01 5 


= h(20) 


h(10) 


= 


0.086 


= h(14) 


h( 


5) 


= 


0.026 


= h(19) 


h(11) 




0.092 


= h(13) 


h( 


6) 


= 


0.038 


= h(18) 


h(12) 


= 


0.094 


= h(12) 



und die Transferfunktion Abb. 7, die den oben erwähnten 
Nachteil des Parzen-Filters nicht hat und deshalb auch für 
Reihen mit Monatswerten verv/enbar ist, wie aus Abb. 8 deut- 
lich wird. 

1.2.3 Ein Vergleich mit einfachen gleitenden 
Durchschnitten 



In der elementaren Zeitreihenanalyse werden zur Bestimmung 
der glatten Komponenten einer Zeitreihe häufig einfache glei- 
tende Durchschnitte verwendet. Die Filterkoeffizienten eines 
solchen gleitenden Durchschnittsfilters sind 

h(k) = 1 , |k| < ^ 

Für die Simulationsreihe ergibt sich Abb. 9 mit einem Theil- 
schen üngleichheitskoef f izienten von 0,03547, der also we- 
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sentlich größer ist als die entsprechenden Koeffizienten der 
"Fenster "-Filter und damit eine merklich schlechtere Schät- 
zung der glatten Komponente als mit jenen Filtern anzeigt. In 
der Tat ermittelt dieser Filter die glatte Komponente als zu 
"glatt". Trotz dieser seit langem bekannten Eigenschaft dieses 
Filters, werden auch heute noch einfache gleitende Durchschnit- 
te z. B. bei Saisonbereinigungsverfahren eingesetzt, obwohl 
dadurch zyklische Prozesse konjunktureller Provenienz entwe- 
der völlig ausgefiltert oder ihre Amplitude zumindest so 
stark abgeschliffen werden, daß sie praktisch nicht mehr in 
Erscheinung treten. Das ist insbesondere dann der Fall, wenn 
solche Durchschnitte mehrfach angewendet werden, wie etwa 
beim weit verbreiteten Census X-1 1 -Verfahren . Es ist in diesem 
Zusammenhang überhaupt kritisch zu fragen, ob die Vorstellung, 
daß die dem Niederfrequenzbereich zuzuordnenden Komponenten 
einer ökonomischen Zeitreihe "glatt" verlaufen, also durch 
Polynome niedriger Ordnung approximierbar seien, historisch 
ihre Wurzel nicht im praktisch ausschließlichen Gebrauch von 
einfachen gleitenden Durchschnitten hat. Diese Vorstellung aber 
ist bis heute eine wichtige Grundlage für die meisten prakti- 
zierenden Saisonbereinigunqsverf ahren . Auf diese interessanten 
Zusammenhänge kann jedoch erst in 3.3. eingegangen werden. 

Der Grund dafür, daß einfache gleitende Durchschnitte im allge- 
meinen "schlechte" Tiefpaßfilter sind, läßt sich leicht an ih- 
rer Transferfunktion erkennen. Für einen symmetrischen gleiten- 
den Durchschnitt lautet diese: 



T 



GD 



m _ sin (XN/2) 

^ ^ N sin(A/2) 



I A I <TT 



Daraus folgt 



lim T 

N-»oo 



GD 



(A) 



r O für A 4 0 
^ 1 für A = O 



KOOPMANS , 19 ) d. h. diese Transferfunktion konvergiert 

punktweise gegen die diskrete Impuls-Antwortfunktion. Diese 
stellt die ideale Transferfunktion eines Filters dar für Korn- 
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ponenten einer Zeitreihe, deren spektrale Masse auschließlich 
in der Frequenz A = O konzentriert ist. Dies ist z. B. dann der 
Fall, wenn eine Zeitreihe nur einen Trend auf weist, aber sonst 
keinerlei niederfrequenten oszillatorischen Komponenten. Hier- 
aus folgt, daß einfache gleitende Durchschnitte nur in diesem 
Spezialfall als adäquate Tiefpaßfilter anzusehen sind. Da sich 
in ökonomischen Zeitreihen in der Regel aber auch der Einfluß 
konjunktureller Prozesse widerspiegelt, ist davon auszugehen, 
daß dieser Spezialfall wohl relativ selten auf tritt, z. B. 
etwa bei Reihen mit Jahreswerten, die zu kurz sind um mehre- 
re Konjunkturzyklen zu umfassen. Anders verhalten sich dage- 
gen die Transferfunktionen derjenigen Filter, die einerseits 
auf dem Optimalprinzip eines minimalen "mittleren quadrati- 
schen Fehlers" beruhen und andererseits durch Normierung und 

Einführung von Fensterfunktionen ( "windowing") ihre endgül- 

1 ) 

tige Gestalt gewinnen. Sie konvergieren mit zunehmender 
Filterlänge N gegen die intendierte Transferfunktion. Prak- 
tisch ist bei diesem Vergleich allerdings eine Einschränkung 
zu machen für kleine Filterlängen, da sich dabei die Schätz- 
ergebnisse bei Verwendung von einfachen gleitenden Durchschnit- 
ten von denen der Fensterfilter numerisch u. U. nur unwesent- 
lich unterscheiden. Die Situation ist vergleichbar mit dem 
Konzept einer konsistenten Schätzfunktion in der Statistik. 
Bekanntlich ergeben sich hierbei im Vergleich zu inkonsisten- 
ten Schätzfunktionen bei kleinen Stichprobenumfängen meistens 
numerisch keine wesentlichen Unterschiede. Trotzdem verzich- 
tet man bei Schätzfunktionen aus guten Gründen nicht auf die 
asymptotische Eigenschaft der Konsistenz. 



1) Da {— !— e^^^,k = O, +1 f±2 , . . . } eine Hilbertbasis im dar- 
/2tt 

stellt, liegt für die Transferfunktionen der Rechteckfilter 
für N oo Konvergenz im quadratischen Mittel gegen die 
ideale Transferfunktion vor. Das gilt auch für die anderen 
Filter, da die zugehörigen Fensterfunktionen für N 
punktweise gegen das Rechteckfenster konvergieren. 




Es sei noch bemerkt, daß sich durch "windowing" nicht nur 
Tiefpaßfilter, sondern auch Bandpaßfilter, Bandsperren usw. 
sowie andere Selektionsfilter in der Regel genau so einfach 
konstruieren lassen, was offensichtlich mit einfachen gleiten- 
den Durchschnitten nicht möglich ist. Voraussetzung ist nur, 
daß sich das einschlägige Integral zur Bestimmung der Filter- 
Koeffizienten geschlossen lösen läßt. 
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1.3. Optimale FIR-Filter 

Außer der oben besprochenen Konstruktion von FIR-Filtern mit 
Hilfe von Fensterfunktionen, sind noch andere Konstruktions- 
prinzipien denkbar, die hier kurz angedeutet werden sollen. 
Kurz deswegen, weil es sehr zweifelhaft erscheint, ob so kon- 
struierte Filter von praktischem Wert für die Analyse ökono- 
mischer Zeitreihen sind. 

Für die Frequenz-Antwortfunktion eines FIR-Filters mit Null- 
phase : 

M . , 

T(X) = Z h(k)e"^ ^ 
k=-M 

wobei die Dauer der Impulsantwortfunktion 2M+1 Zeiteinhei- 
ten beträgt, kann geschrieben v/erden: 

M 

T (A) = h (o) + 2 Z h (k) cos Ak 
k=1 



da gilt 

h(k) = h(-k) 

Es soll nun ein Tiefpaßfilter konstruiert werden mit folgen- 
der Spezifikation: 



iT(A) 
iliniimn uii 












imrrrTTTTTrm 

I 

I 

f 1 

1 


! i n 1 iiiiuiii lu.iin / 






^2 mini ! n ////// /////// //' 





d. h. im Paßband O ^ A ^ A^ soll die ideale Transferfunk- 
tion (die dort den Wert Eins annimmt) mit einem maximalen 
Fehler von 6^ approximiert werden und im Stoppband 
A^ £ A £ TT mit einem maximalen Fehler von 6^ (idealer Wert: 
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Null). Das Übergangsband (A^, A^) bleibt unspezif iziert 

Die Parameter M, 6 ^, 62 ^ A^ und A 2 können nicht unabhängig 

voneinander festgelegt werden. Beim Konstruktionsverfahren 
von (HERMANN-SCHUESSLER, 14) können M , 6 ^ und 62 festge- 
legt werden, eine a pr ior i-Kontrolle über die praktisch oft 
wichtigen Parameter A^ und A 2 ist also nicht möglich. Dage- 
gen ist z. B. bei dem Verfahren von (PARKS-McCLELL AN , 27) 

die Möglichkeit gegeben, neben M auch A^ und A 2 festzulegen. 
Für die maximalen Fehler 6 ^ und 6 ^ sind dann allerdings im 
Allgemeinen keine a priori-Restr ikt ionen zu realisieren. 

Bei dem eben erwähnten Verfahren von HERMANN- SC HUESSLER 
wird eine Approximation an eine vorgegebene ideale Trans- 
ferfunktion versucht, die als "equiripple approximat ion" 
bekannt ist. Die Anzahl der lokalen Maxima und Minima der 
Approximationsfunktion im Frequenzbereich 0 A £ tt hängt 
dabei von M, also der Dauer der Impulsantwortfunktion ab, 
denn für T(A) kann geschrieben werden: 

^ k 

T ( A ) = 1 Uj, (cos A ) 

k=o 

wobei Konstanten sind, die von h(k) abhängen. Somit ist 
T(A) ein trigonometrisches Polynom der Ordnung M, das im 
Intervall O < A < tt höchstens M-1 lokale Extremv/erte besitzt. 
Nimmt man das abgeschlossene Intervall [0,tt], dann kann 
gezeigt werden, daß T(A) bei A = 0 und A = tt entweder ein 
xMaximum oder ein Minimum annimmt. Somit v/eist T ( A) im ab- 
geschlossenen Intervall höchstens M+1 lokale Extrema auf. 
Liegen Ep Extrema im Passband und E^ Extrema im Stoppband, 
dann gilt: 

Ep + Eg = M + 1 

Es läßt sich nun ein nicht-lineares Gleichungssystem mit 
2M Gleichungen angeben zur Bestimmung der unbekannten M+1 



1) Für die Ausführungen in diesem Abschnitt vgl. (OPPENHEIM- 
SCHÄFER, 25 , S. 255 ff) 
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Filterkoeffizienten und der M-1 Frequenzen, bei welchen die 
M+1 Extrema liegen (bei A = 0 und A = tt liegen immer Extrema) 

Zweifellos liegt der größte Nachteil des eben in seinen Grund 
Zügen skizzierten Verfahrens darin, daß die Größen A^ und A^f 
die im praktischen Fall nach substanzwissenschaftlichen Über- 
legungen festgesetzt werden müssen, nur zufällig realisiert 
werden können. 

Dieser Nachteil tritt beim Verfahren von PARKS-McCLELLAN 
nicht auf. Approximationstheoretisch stellt es eine Tscheby- 
scheff -Approximation dar (vgl. dazu ACHIESER, 1 ). Der Ap- 

proximationsfehler ist dabei folgendermaßen definiert: 

F (A) : = [T* (A) - T(A) ] G(A) 

Dabei ist T*(A) die vorgegebene (ideale) Transferfunktion, 
T(A) die realisierte Transferfunktion und G(A) eine vorge- 
gebene Gewichtsfunktion, die es gestattet, Pass- bzw. Stopp- 
bänder in ihrer Bedeutung für das zu untersuchende Problem 
unterschiedlich zu akzentuieren. Das Verfahren besteht nun in 
einem Algorithmus, der den Ausdruck 

max I F ( A) I 
0 £ A < A^ 

minimiert. Dieser beruht wesentlich auf dem sog. "Alternan- 
tensatz". Für Einzelheiten sei auf die zitierte Originalar- 
beit verwiesen. 

Praktische Versuche mit so konstruierten Filtern führten je- 
doch zu enttäuschenden Resultaten, trotz der Möglichkeit, 
über G(A) Einfluß auf das Resultat der Approximation zu neh- 
men. Versucht wurde z. B. die Konstruktion eines Tiefpaßfil- 
ters mit dem Passband O £ f £ f ^ = 0,05 und dem Stoppband 
0.053 = ^2 ^ ^ Filterlängen 23, 51 , 99. Dabei 

wurden Passband und Stoppband in der Relation 5 : 1 bzw. 

10 : 1 gewichtet. Wie die Abb. 10-15 zeigen, verlaufen die 
erzielten Transferfunktionen völlig unbefriedigend. Versucht 
man etwa damit, eine konkrete Reihe zu filtern, dann er- 




23 



geben sich ebenfalls völlig indiskutable Resultate. Die Abb. 
16-18 zeigen die Anwendung dieser Filter auf die Reihe "Netto- 
produktion der verarbeitenden Industrie". Zunächst fällt auf, 
daß die gefilterten Reihen systematisch über oder unter der 
Originalreihe liegen. Das liegt daran, daß die entsprechen- 
den Transferfunktionen für A = 0 entweder größer als Eins oder 
kleiner als Eins sind. Da ökonomische Reihen meistens einen 
Trend aufweisen, führt jede Abweichung |T(o) - 1 |, auch wenn 
sie numerisch klein ist, zu systematischen Verfälschungen. 

Da es aber unmöglich ist, a priori sicherzustellen, daß 
T (0) = 1 ist, scheinen die Einsatzmöglichkeiten derartiger 
Filter als Tiefpaßfilter für ökonomische Zeitreihen recht ge- 
ring zu sein. Aber selbst wenn für eine konkrete Problemstel- 
lung dieser Defekt unwesentlich sein sollte, ist an eine prak- 
tische Verv/ertbarkeit so konstruierter Filter kaum zu denken, 
da sie schlechte Ausblendeigenschaften aufweisen, v/ie aus den 
abgebildeten Transferfunktionen leicht ersichtlich ist. Dies 
gilt selbst für den Fall großer Fi It er längen , wie sie in praxi 
nie realisiert werden können. 
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2. Rekursive Filter 



2.1. Einfache rekursive Filter 



Die einfachsten rekursiven Filter erhält man aus der allge- 
meinen Filtergleichung 

M N 

Z a, y(t-k) = ^ b x(t-r) , a = 1 

k=0 r=0 ^ 

wenn M=1 und N=0 gesetzt wird.. Dann ergibt sich: 

y(t) = b^ x(t) - a^ y(t-1) 

Für geeignet gewählte Koeffizienten erhält man daraus leicht 

1 ) 

z. B. einen Tiefpaßfilter . Auch läßt sich unschwer sicher- 
stellen, daß die zugehörige Transf erf unktion bei A = 0 genau 
den Wert Eins annimmt. Es ist: 



T (z) 



N 

E b^ 
r=0 
M 

Z 

k=0 



b^ z 

0 • 



^k ^ 



-k 



1 - a^ z 



-1 



Aus der Bedingung T(1) = 1 ergibt sich dann 

Dieser Filter hängt somit nur von einem Parameter ab. Für po- 
sitives a^ liegt ein Tiefpaßfilter vor, für negatives a^ ein 
Hochpaßfilter. Dabei verlaufen die Transferfunktionen umso 
schärfer, ja näher |a^| bei Eins liegt (aus Stabilitätsgrün- 
den ist |a^ I = 1 auszuschließen) . Die Abb. 19 zeigt die zum 
"half -pov/er "-Punkt f^=0,05 gehörige Transferfunktion. 

So konstruierte Filter sind natürlich sehr einfach zu hand- 
haben, haben allerdings den Nachteil, daß ihre Transf er funk- 
tionen relativ unscharf sind, obv7ohl sie unvergleichlich 
schärfer verlaufen als etwa diejenigen der unter großem 



1) Für nachfolgende Ausführungen vgl. (CADZOW, 6). 
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theoretischen und praktischen Aufv/and erzielten FIR-Filter. Es 
ist insbesondere zu beachten, daß sie auch bei kurzen Zeit- 
reihen einsetzbar sind. Trotzdem werden sie sich je nach Pro- 
blemstellung evtl, als ungeeignet erweisen, z. B. wenn eine 
Zeitreihe außer einem Trend noch eine nennenswerte konjunk- 
turelle Komponente besitzt. Es besteht dann nämlich die Ge- 
fahr, daß diese unter Umständen sehr stark gedämpft v/ird , 
ein Problem, das schon bei den einfachen gleitenden Durch- 
schnitten auf trat. 

Verschärfungen einer Transferfunktion können aber erzielt 
werden, wenn der zugehörige Filter in Serie geschaltet wird. 
Das eben erwähnte Problem wird dann in einer Relevanz redu- 
ziert, da sich die Durchlass- und Ausblendeigenschaften mit 
zunehmender Filterordnung verbessern. Wird etwa der obige 
Tiefpaßfilter iterativ n-mal angewandt, dann lautet seine 
Transf er funkt ion 

(1-a ) z " 

T„(z) = ^ 

n z-a^ 

die in Abb. 20 für n = 2,45 v/iedergegeben ist (f^=0,05). 

Je größer allerdings n ist, umso größer ist natürlich die zu- 
gehörige Phase, Was unter Umständen unervainscht sein kann. 

Die erwähnten Effekte kommen in den Abb. 21-23 zum Ausdruck, 
wo die schon in 1.1.1. erwähnte Simulationsreihe durch einen 
Tiefpaßtf ilter geschickt wurde mit n = 1, 2, 45. 

Eine merkliche Phasenverschiebung ist hier jedoch erst bei 
großem n festzustellen. Offensichtlich liefert der Filter hier 
auch bei kleinem n recht zufriedenstellende Resultate. 

Für viele Probleme der Zeitreihenanalyse, z. B. der Isolation 
der konjunkturellen Komponente einer Zeitreihe, sind Bandfil- 
ter von Bedeutung. Ein einfacher Bandfilter wird z. B. durch 
folgende Transferfunktion gegeben: 
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T(z) = 



ßz" 

(z - re^\) ) (z - re 



ßz" 



2 2 
z - (2r cos Aq ) z+r 



Dabei ist 

ß = (1 - r) |e^^o - re | 



ein Normierungsfaktor, der sicher stellt , daß T(A) = 1 ist für 

A = Ao . 

Die Bandbreite dieses Filters ergibt sich dabei aus der Be- 
dingung 



|t(A) 




was zu einer quadratischen Gleichung in cos A führt: 
2 

2 ( 1 +r ) cos A 

cos A - [ -] cos A 



+ 



^ ^ — [4r cos 2 A 
4r^ o 



( 1 -r) ^ ] + cos^ A = 0 
o 



mit den beiden Lösungen: 



cos A^ 
cos A^ 



(1+r^) cos A^ (1-r) 

2r “2r 

(1+r^) cos A^ (1-r) 

2r 2r 



+r ) ^- ( 1 -6r+r^) cos^A^ 
V(1+r) (l-6r+r^)cos^A 

o 



woraus sich A^ und A^ ergeben. 

Praktisch wird allerdings üblicherweise umgekehrt vorgegan- 
gen: Man gibt die "half-power "-Punkte A^ und A 2 vor (defi- 
niert also dadurch die gewünschte Bandbreite) und bestimmt 
damit r und A^, wobei sich r als diejenige Wurzel des Poly- 
noms 

+ (a-ß)r^ - 8ar^ + (14a-2ß-2) r"^ 



- 8ar + (a-ß) 



r^ ^ 1 



O 
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ergibt, die im Intervall (0,1) liegt mit 

2 

a = (cos + cos X^) 

2 

B = (cos X^ - cos X^) 

X^ folgt dann aus der Beziehung 

cos A = — ö— (cos A. + cos A,.) 

° 1 + 1 2 

Abb.24 zeigt die Transf erf unktion dieses Filters für f^=0,02 
und f2=0,04. Verschärfungen sind wiederum zu erzielen durch 
Serienschaltung, v/ie in der Abb. 25 dargestellt ist. 

Versucht man z. B. in obiger Simulationsreihe das niederfre- 
quente Signal 200 sin 0.06 irt mit diesem Filter auszufil- 
tern, dann zeigen die Abb. 26-28, daß es dem Filter bei klei- 
ner Ordnung (und damit auch mit geringer Phasenverschiebung) 
nur unzureichend gelingt, eine Trennung von Trend und zykli- 
scher Komponente herbeizuführen. In etwa gelingt das erst z. 

B. bei n = 4 und relativ engem Band [0.025, 0.035], aller- 
dings mit erheblicher Phasenverschiebung. Man könnte nun ver- 
sucht sein, die unzureichende Trennschärfe zwischen diesen bei- 
den Komponenten dadurch zu erhöhen, daß man eine sehr kleine 
Bandbreite wählt und eine mäßig hohe Filterordnung. Wie Abb. 31 
allerdings zeigt, führt ein derartiges Vorgehen jedoch schnell 
zu numerisch unsinnigen Ergebnissen. Es erscheint höchst 
zweifelhaft, ob es gelingt, das alte Problem einer Trennung 
von Trend- und Konjunkturkomponenten mit Hilfe eines derar- 
tigen einfachen Bandfilters zu lösen. Für Zwecke einer Konjunk- 
turdiagnose und Konjunkturprognose wäre jedoch eine solche Tren- 
nung auf filtertheoretischer Basis höchst wünschenswert, um 
damit endlich die bisherigen ad-hoc-Methoden abzulösen (v7ie 
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z. B. Berechnung von Wachstumsraten, gleitende Regressionen 
usw. ) 



2.2, Optimale rekursive Filter 
Genau v;ie bei den FIR-Filtern lassen sich auch bei den II R- 
Filtern Konstruktionsprinzipien angeben, die es gestatten, 
"optimale" Filter in dem Sinne zu konstruieren, daß ein vor- 
gegebenes Fehlerkriterium erfüllt v/ird. Dabei wird so vorge- 
gangen, daß die gev.ünschte Transfer funktion eines Filters an 
den Stellen X^, i = 1,2,... M vorgegeben wird. Wählt man 
V7ie z. B. (STEIGLITZ, 36' den mittleren quadratischen Fehler 
als Kriterium, dann ergibt sich für diesen der Ausdruck: 

M - ? 

MQF = E [ 1t(X ) | - |T(A )p 
i=1 ^ ^ 



mit 



T(X) : = realisierte Transferfunktion 

T(A) : = erwünschte Transfer funktion 

wobei T(X) zweckmäßigerweise in Serienschaltung von Filternder 
Ordnung 2 realisiert v;ird: 



T (z) 



K 

A n 

k=1 



1 +a z 
k 






-2 

z 






K ist dabei die Ordnung des Filters. MQF ist eine Funktion der 
4K Parameter (a^ , b^, c^, d^, ... a^^, b^^, c^^, dj^) . Die Mini- 

mierung von MQF führt zu einem System von 4K nicht-linearen 
Gleichungen, das z. B. mit der Methode von FLETCHER-POWELL 
gelöst werden kann. 



1) Es soll hier nicht unerwähnt bleiben, daß dies in der Zwi- 
schenzeit mit Hilfe einer speziellen Kombination von rekur- 
siven Filtern gelungen ist, die aber hier nicht dargestellt 
werden soll, da ihre Entwicklung erst nach Abschluß des 
vorliegenden Projektes erfolgte. Es sei deshalb auf weite- 
re Publikationen des Autors verwiesen. 
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Versuche, mit so konstruierten reluir^iven Filtern zu arbeiten, 
führten zu Transf erf unktionen , die auch bei relativ niedriger 
Ordnung als für viele praktische Zwecke hinreichend scharf be- 
zeichnet werden können, ganz im Gegensatz zu den optimalen 
FIR-Filtern. Allerdings sind auch sie nicht als problemlos an- 
zusehen, da genau wie bei den optimalen FIR-Filtern keine Kon- 
trolle über den Verlauf einer realisierten Transferfunktion bei 
A = 0 möglich ist, wenn auch die Differenz |t( 0) - 1 | prak- 

tisch immer sehr viel geringer ist als bei den FIR-Filtern. 
Hinzu kommt noch, daß sich je nach Filterordnung unter Umstän- 
den erhebliche Phasenverschiebungen bemerkbar machen, die bei 
gewissen Problemstellungen sicherlich störend in Erscheinung 
treten. Generell könnte man diese optimalen rekursiven Filter 
als bedingt tauglich für die Analyse ökonomischer Zeitreihen 
ansehen. Da auf solche Filter im Zusammenhang mit der Dar- 
stellung eines Bereinigungsverfahrens in Abschnitt 3.5 zurück- 
zukommen sein v/ird , soll hier auf die Wiedergabe von Beispie- 
len verzichtet werden. 

2.3. Der Kerbenfilter 



Neben den einfachen rekursiven und den optimalen rekursiven 
Filtern gibt es noch solche, die auf "direktem" Weg konstru- 
iert werden können. Mit "direkt" ist dabei gemeint, daß der 
Filter durch Festlegung der Null-Stellen und Pole der Trans- 
ferfunktion realisiert wird. Schreibt man die Transf er funk- 
tion eines rekursiven Filters in der Form: 

b (z - z ) (z - z ) ... (z - z ) 

T(z) = ' 1 

(z - ) (z - p^) ... (z - p^) 

wobei z^ (i = 1 , 2, . . . M) die Nullstellen und 

Pj ( j = 1 / 2, . . . N) die Pole von T(z) sind, dann ist leicht 

einzusehen, daß sich durch eine geschickte Wahl der z und 
Pj eine erv,ünschte Transferfunktion evtl, konstruieren läßt. 
Dabei hat b den Charakter einer Normierungskonstant.en , die 
so gewählt werden kann, daß T (z ) an einer vorgegebenen Stel- 
le einen vorgegebenen Wert erreicht. 
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Für die Lösung praktischer Probleme der Zeitreihenanalyse 
hat sich der Einsatz des sog. "Kerbenfilters" als besonders 
zweckmäßig herausgestellt. Darunter ist ein Filter zu ver- 
stehen, dessen Transferfunktion idealerweise folgende Ge- 
stalt hat: 




o 



Eine solche Transferfunktion läßt sich nun mit hoher Ge- 
nauigkeit durch folgende Überlegungen gev/innen: 

Da bei ^ einzige Nullstelle liegen soll, muß auf 

alle Fälle gelten: 

iA o 

z^ = p 

Wählt man nun den ersten Pol in der Nähe dieser Nullstelle, 
also , dann ist der Quotient 



Z - Pi 

für alle z 4 e^^^ außerhalb des "Kerbenbereiches" liegen 

etwa gleich Eins. Eine solche Transferfunktion der nieder- • 
sten Ordnung (M = 1 , N = 1 ) führt aber zu komplexen Filter- 
koeffizienten, die praktisch nicht brauchbar sind. Reelle Ko- 
effizienten erhält man erst, wenn man mindestens die Ordnung 
(M = 2, N = 2) wählt, indem man setzt: 



z 



2 



z 



1 



P2 Pi 
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Die Transfer funktion lautet dann: 



T(z) 



, iXOv / -iAox 

( z - e ) (z - e ) 

(z - ) ( Z - ) 



z 



iA 

e 



Diese Transferfunktion verläuft allerdings noch ziemlich unbe- 
friedigend, insbesondere ist die Kerbe nicht hinreichend symme- 
trisch um Aq. Eine Verbesserung der Verlaufsform kann nun v/ie- 
der wie bei allen Filtern durch eine Serienschaltung erzielt 
werden. Allerdings empfiehlt es sich dabei, eine geringe Pol- 
variation bei den sukzessiven Filtern vorzunehmen, wobei die 
Nullstelle aber jeweils gleich bleibt. Eine solche Filterfol- 
ge wurde von (WEAVER et al., 39) für Elektrokardiagramm-Ana- 
lysen benutzt. Dabei wurde die Filterordnung (M = 6, N = 6) 
gewählt und für die Pole eine "symmetrische Anordnung" vorge- 
schlagen. Die Pol-Nullstellen-Konstellation ist nachstehend 
wiedergegeben für 0 £ ^H/2: 




Für die Pole gilt daher: 



= cos Aq - A cos (Aq - a) + i [sin A^ - A sin (A^ -a) ] 

P2 = (1 - A) cos Aq + i (1 - A) sin A^ 

p^ = cos A^ - A cos (A^ + a) + i [sin A^ - A sin (A^ + a) ] 

(Vgl. dazu CADZOVJ, 6 ) . 
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Die Pole hängen von den Größen a und A ab. Für die Transfer- 
funktion des Kerbenfilters ergibt sich dann insgesamt der 
Ausdruck: 



, ^iAOv3 / -iAo^3 

T (z ) = ^ ) (z - e )__ 



3 

n ( z 

i=1 



Z = e' 



iX 



p^) ( z - p^) 



Die Normierungskonstante wird zweckmäßigerweise so festge- 
legt, daß T(A) den Wert Eins annimmt an der Stelle A = 0, 
d. h. also: ^ 

n (1 - p^) (1 - p.) 

^ ^ i^J 

iAO\3 f ^ — iAO\3 

Ist a festgelegt, dann hängt die Transferfunktion nur noch von 
A ab. Je kleiner A ist, umso enger ist die Kerbenöffnung. 
(WEAVER et al., 39) setzen a = 72^. Die Abb. 32,33 zeigen die 
Transferfunktion des Kerbenfilters für A = 0.01 und A = 0.1 
und Aq = 2tt/1 2 . 

Wie diese Abbildungen zeigen, entsprechen die so gefundenen 
Transferfunktionen noch keinesfalls den theoretischen Vor- 
stellungen, da sie in der Nähe der Kerbenfrequenz A^ "Nasen" 
aufweisen. Es bleibt deshalb die Aufgabe, diesen Defekt zu be- 
seitigen durch eine Veränderung der Lage der Pole. Darauf wird 
noch zurückzukommen sein. 
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3. Die Verwendung von rekursiven 
Filtern zur Lösung praktischer 
Probleme der Zeitreihenanalyse 

Nachfolgend werden nun einige Anwendungsmöglichkeiten rekursi- 
ver Filter dargestellt, wobei aber die im Projekt entwickel- 
ten Verfahren nur skizziert werden, da diese in speziellen Ver- 
öffentlichungen ausführlich erörtert werden. Außerdem waren 
diese Verfahren bei Abschluß des Projektes noch nicht ganz an- 
wendungsreif bzw. ausgetestet. 

Im einzelnen wird zunächst über ein Verfahren berichtet, das 
als geeignetes Analyseinstrument anzusehen ist, um die alte 
Frage nach der Existenz bzw. Nichtexistenz langer Konjunktur- 
wellen definitiv zu beantworten. Im Anschluß daran wird auf 
das vieldiskutierte Problem der Saisonbereinigung ökonomischer 
Zeitreihen eingegangen. Dazu werden zu Beginn grundsätzliche 
Bemerkungen gemacht, um dann in einer längeren Untersuchung 
nachzuweisen, daß eine objektive Verknüpfung von "Zeitbereichs- 
analyse" und "Frequenzbereichsanalyse" nicht möglich erscheint. 
Daran anschließend wird gezeigt, warum eine Verwendung von re- 
kursiven Filtern gerade für Zwecke der Saisonbereinigung sinn- 
voll ist. Schließlich werden zwei neue Verfahren skizziert. 

Auf die Wiedergabe von Ergebnissen wird dabei allerdings für 
Verfahren I verzichtet, da ihre Diskussion umfangreiche Ablei- 
tungen und Überlegungen erfordert, auf die hier aus Platzgrün- 
den verzichtet werden muß. Hier muß wieder auf die schon ein- 
gangs erwähnten, sich in Vorbereitung befindlichen, Veröffent- 
lichungen hingewiesen werden. 

3.1. Zum Problem der "long Swings " 

Die Diskussion über die Existenz bzw. Nichtexistenz langer oder 
längerer konjunktureller Wellenbewegungen hat in jüngster Zeit 
wieder an Intensität gewonnen, aus Gründen, die hier nicht dis- 
kutiert werden sollen. Ausgangspunkt ist folgender Sachverhalt: 
Wirtschaftsforscher und -historiker unterschiedlichster Prove- 
nienz formulierten die Hypothese, daß sich ökonomische Ent- 
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Wicklungsprozesse nicht "geradlinig", sondern in Wellenbewe- 
gungen vollzögen, deren Zyklusdauer aber von den einzelnen 
Forschern aufgrund ihrer empirischen Untersuchungen durchaus 
verschieden veranschlagt wird. Zw^ingt schon dieser Sachver- 
halt zu grundsätzlicher Skepsis gegenüber den vorgelegten 
Resultaten, so ist darüber hinaus ganz generell zu fragen, 
ob denn das adäquate Instrumentarium vorhanden ist, um der- 
artige Hypothesen überhaupt prüfen zu können. 

Bemüht man Methoden der traditionallen Zeitreihenanalyse, 
dann wird davon ausgegangen, daß sich eine Reihe a priori 
in eine gewisse Anzahl von Komponenten zerlegen lasse, un- 
ter denen sich möglicherweise auch eine "zyklische" befin- 
de, die aber, obgleich sie im konkreten Fall durchaus mit 
bemerkenswerter Regelmäßigkeit verlaufen könne, optisch 
meistens nicht oder nur sehr schwach in Erscheinung trete, 
da andere Komponenten dominierten, insbesondere der "Trend". 
Gelingt es nun, den Trend zu eliminieren, dann kann prinzi- 
piell eine vorhandene zyklische Komponente isoliert darge- 
stellt werden. Aber v/ie soll der Trend eliminiert werden? 

Die traditionelle Analyse bedient sich dazu verschiedener 
Methoden. Einmal wird von der Vorstellung ausgegangen, die 
Trendkomponente ließe sich durch ein Polynom dar stellen, 
dessen Koeffizienten geschätzt werden können, v/as dann zu 
einer verhältnismäßig einfachen "Ausschaltung" des Trends 
führt. Da aber der dabei vorauszusetzende Polynomgrad sub- 
stanzwissenschaftlich praktisch nie bestimmbar ist, ergeben 
sich Resultate, die nicht frei von Willkür sind. ANDERSON 
sen. hat das anhand der Untersuchungen von KONDRATIEFF, der 
in der geschilderten Weise vorgegangen ist, mit aller Deut- 
lichkeit gezeigt (vgl. ANDERSON sen., 2 ) . Es bedarf also 
keiner Diskussion, daß eine derartige Vorgehensweise als me- 
thodisch inadäquat zurückzuweisen ist. 
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Zum anderen verwendet die traditionelle Analyse aber auch Ver- 
fahren, die in mancher Hinsicht als wesentlich "flexibler" als 
die Vorgabe eines "starren" Trendpolynoms anzusehen sind. Hier 
unter fallen alle Prozeduren, die unter dem Oberbegriff "glei- 
tende Durchschnitte" subsumiert werden können. Spätestens seit 
(SLUTZKY, 35) ist jedoch bekannt, daß solche Durchschnitte 
unter Umständen zu erheblichen Fehlschlüssen verleiten können, 
da sie möglicherweise Zyklen induzieren, die im originären 
Datenmaterial, d. h. in einer ungeglätteten Zeitreihe, nicht 
enthalten sind. Als Paradebeispiel sei hier auf KUZNETS 
verwiesen. Mit den Hilfsmitteln der modernen Zeitreihenana- 
lyse kann nämlich gezeigt werden, daß die von KUZNETS heraus- 
gearbeitete Zyklendauer das Resultat seiner Glättungsprozedu- 
ren ist (vgl. FISffi^AN, 8 ) . Somit sind also auch derartige 
Methoden der vorliegenden Problemstellung nicht angemessen. 

Die moderne Zeitreihenanalyse - und hier sei insbesondere 
an die Spektralanalyse gedacht - vermeidet nun eine a priori- 
Komponentenzer legung einer Zeitreihe. Durch eine Analyse, 
die prinzipiell im Frequenzbereich durchgeführt wird, ergibt 
sich die Möglichkeit einer Zeitreihenzerlegung, ohne daß eine 
vorherige Elimination gewisser Komponenten erforderlich ist. 
Leider ist dieser unbestreitbare Vorteil einer Frequenzbe- 
reichsanalyse an eine entscheidende Voraussetzung gebunden, 
nämlich der (schwachen) Stationär ität des zugrundeliegenden 
stochastischen Prozesses. Dies impliziert vor allem, daß eine 
Zeitreihe keine evolutor ischen Bewegungen aufweist, was aber 
gerade typisch für ökonomische Reihen ist. Dies führt dazu, 
daß infolge von " leakage"-Ef f ekten das Auflösungsvermögen 
von Spektralfenstern im Niederfrequenzbereich (und dieser 
interessiert ja nur für das vorliegende Problem) so gering ist 
daß ein geschätztes Spektrum ohne jegliche Information bezüg- 
lich der Existenz von längeren zyklischen Bewegungsabläufen 
ist. Will man informative Spektren, dann ist es unumgänglich. 




36 



vor der eigentlichen Spektralanalyse eine Trendbereinigung 
durchzuführen. Genau v;ie bei den traditionellen Verfahren 
stellt sich auch hier sofort die Frage nach dem "Wie?". 

Üblicherweise wird genauso vorgegangen wie bei jenen, d. h. 
also Polynomansatz und/oder gleitende Durchschnitte. Aber 
auch hier kann ein Polynomgrad nicht objektiv festgelegt wer- 
den, so daß insgesamt die Situation nicht als günstiger zu 
beurteilen ist als bei den herkömmlichen Verfahren. Grundsätz- 
lich vorteilhafter scheint sie allerdings bei Verwendung von 
gleitenden Durchschnitten zu sein: Solche Filter können be- 
züglich ihrer Effekte anhand der Transferfunktion im Frequenz- 
bereich sehr detailliert untersucht werden, V 7 as mit den Werk- 
zeugen der traditionellen Analyse im allgemeinen nicht möglich 
ist. Es stellt sich deshalb die Frage, ob Filter existieren, 
die so beschaffen sind, daß eine Trendelimination möglich ist, 
ohne daß dadurch unerwünschte Nebeneffekte auf treten. Erst 
dann v/äre eine informative Spektralanalyse möglich. 

Wie muß nun ein Trendfilter beschaffen sein, damit Zeitreihen 
mit Hilfe der Spektralanalyse im Hinblick auf die Existenz 
niederfrequenter Zyklen adäquat untersucht werden können? Die 
Antwort lautet: ein solcher Filter muß so beschaffen sein, daß 
seine Transferfunktion T(X) zv/ei Bedingungen erfüllt. Diese 
sind (vgl. HOWREY, 15): 

1) T (Aq) =1 

2) T' (Aq) = 0, Aq€(0,tt) 

wenn die Hypothese getestet werden soll, daß eine Reihe eine 
zyklische Komponente der Frequenz A^ enthält . Durch die Er- 
füllung der ersten Bedingung wird sichergestellt, daß weder 
Verstärkungs- noch Abschwächungseffekte eintreten, d. h. also, 
eine vorhandene bzw. nicht vorhandene Spitze im Spektrum an 
der Stelle A = A^ ist nicht das Resultat der vorgenommenen 
Filteroperation. Ist zudem die zweite Bedingung erfüllt, dann 



1) Selbstverständlich muß sie außerdem im "Nullband" den Wert 
Null annehmen. 
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ist man sicher, daß eine eventuell im Spektrum auf tretende 
Spitze an der "richtigen" Stelle sich befindet, d. h. also, 
daß Verschiebungen im Frequenzbereich und damit ein Ausweis 
von zu großen Periodenlängen (wenn T' (Xq)< 0) oder zu kleinen 
Periodenlängen (v/enn T'(A^)>0) vermieden v/ird. Es ist fast 
überflüssig zu bemerken, daß beide Bedingungen nicht nur für 
den Frequenzpunkt X = erfüllt sein müssen, sondern auch 
für eine genügend breite Umgebung dieses Punktes, da je Zy- 
klen-Hypothesen keine exakten Periodizitäten, sondern nur 
durchschnittliche Periodizitäten implizieren, eine Vorstel- 
lung, die bekanntlich der gesamten Spektralanalyse zugrunde 
liegt . 

Schon vor einigen Jahren haben (KÖNIG und WOLTERS, 20 ) die 
Frage untersucht, welche der bisher üblichen gleitenden 
Durchschnittsfilter die oben formulierten Bedingungen erfül- 
len. Sie kamen zu dem Schluß, daß sie bei keinem erfüllt 
sind. Eine gewisse Ausnahme bildete lediglich ein spezieller 
Hochpaßf ilter^ der allerdings nur zum Test von Zyklen ge- 

eignet ist, deren Periodenlänge kleiner als zehn Jahre sind. 
Die eigentlichen " long-swing" -Hypothesen implizieren aber 
eine wesentlich längere Zyklendauer (z. B. KUZNETS 20 Jahre, 
KONDRATIEFF 50 Jahre etc.), so daß auch dieser Filter eigent- 
lich nicht in Frage kommt. 

Ein für vorliegende Problemstellung geeigneter Trendfilter 
wäre zweifellos der oben diskutierte Kerbenfilter, wenn es 
gelingt, seine "Nasen" zu beseitigen, so daß die erwähnten 
Forderungen erfüllbar werden. Die Lösung dieses Problems so- 
wie detaillierte Ausführungen zur Verwendung dieses Filter- 
typs nebst Anwendungsbeispielen findet sich in einer Arbeit 
von (SCHULTE, 34) , so daß sich weitere Erörterungen hier er- 
übrigen . 
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3.2. Grundsätzliche Probleme der 

Saisonbereinigung ökonomischer Zeitreihen 



3.2.1 . Vorbemerkungen 

Das Problem der Saisonbereinigung ökonomischer Zeitreihen 
wird schon seit gut 1CX) Jahren diskutiert. Diese Diskussion 
hat in den letzten Jahren erheblich an Intensität gewonnen. 

Das hängt nicht nur damit zusammen, daß in diesem Zeitraum 
neue Saisonbereinigungsverfahren entwickelt und auch prak- 
tiziert wurden, sondern auch damit, daß die methodischen 
Grundlagen dieser Verfahren vergleichend kritisch analysiert 
wurden und noch werden. Während sowohl die Möglichkeit von 
sinnvollen Saisonbereinigungen als auch die Notwendigkeit 
derartiger Bereinigungen für die Lösung konkreter Proble- 
me der Wirtschaftsforschung und der Wirtschaftpolitik un- 
strittig sind, gehen die Meinungen darüber, wie die jewei- 
ligen Vor- und Nachteile einzelner Verfahren zu bewerten 
bzw. als solche überhaupt objektiv zu erkennen sind, z. 

Zt. noch im allgemeinen nicht unerheblich auseinander. Ver- 
mutlich ist es nur so zu erklären, daß verschiedene Insti- 
tutionen wie etwa Wirtschaftsforschungsinstitute, das Sta- 
tistische Bundesamt usw. immer noch mit unterschiedlichen 
Verfahren arbeiten und demzufolge u. U. zu erheblich diver- 
gierenden Schlüssen kommen (Vgl. dazu etwa die mit verschie- 
denen Verfahren saisonbereinigte Reihe "Arbeitslose" bei 
KARR, 17). Angesichts der Tatsache, daß manche wirtschafts- 
statistische Zeitreihe in ihrer saisonbereinigten Form heu- 
te von beträchtlicher politischer Brisanz ist, stellt die- 
ser Sachverhalt nicht gerade ein Ruhmesblatt für die Statistik 
dar, nährt es doch u. U. das weitverbreitete Vorurteil, daß 
man mit Statistik schließlich alles beweisen kann. 
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Die vorliegende Abhandlung verfolgt nicht das Ziel, die me- 
thodischen Grundlagen und Konzepte verschiedener Verfahren 
vergleichend zu analysieren, obwohl es sich für manche zu 
diskutierenden Aspekte als unumgänglich erweisen wird, auf 
die Konstruktionsprinzipien konkreter Verfahren einzugehen, 
wenn auch manchmal nur zum Zweck der Illustration der hier 
vorzutragenden Gedankengänge. Auch wird nicht der Versuch ge- 
macht, über numerische Unterschiede zwischen mit alternativen 
Verfahren bereinigten Zeitreihen zu diskutieren, etwa mit dem 
Ziel, herauszufinden, welches Verfahren wohl die "wahre" Sai- 
sonkomponente am "besten" schätzt oder der "wahren" saisonbe- 
reinigten Reihe am "besten" entspricht. Solche Vergleiche, 
die auf den einschlägigen Fachtagungen oft Gegenstand heißer 
Kontroversen sind, erweisen sich nämlich bei näherer Betrach- 
tung als nicht sehr sinnvoll, wie noch zu zeigen sein wird. 
Vielmehr sollen hier nur "Grundprobleme" der Saisonbereini- 
gung behandelt werden, also Probleme, über die vor der Kon- 
struktion eines Verfahrens Klarheit bestehen muß, sollen nicht 
spezielle heute zur Verfügung stehende oder noch zu entwickeln 
de zukünftige Verfahren Anlaß zur Entstehung immer neuer Pro- 
bleme (oder Scheinprobleme?) geben, die zu nicht endenden 
Grundsatzdiskussionen führen. Daran anschließend v/erden Vor- 
schläge zur Konstruktion von Saisonbereinigungsverfahren unter 
breitet, von denen der Verfasser hofft, daß diese nicht nur 
zu einer Vereinheitlichung der gegenwärtigen Diskussion des 
Problemkomplexes "Saisonbereinigung" führen, sondern mögli- 
cherweise auch zu einer Vereinheitlichung der Praxis der Sai- 
sonbereinigung, was wohl angesichts des oben skizzierten Sach- 
verhaltes sicherlich wünschenswert wäre. 




40 



3.2.2. Zum Vergleich von Saisonbereinigungsverfahren 

Die Grundhypothese, die hinter allen bisherigen Saisonberei- 
nigungsverfahren steht, ist die Auffassung, daß sich eine 
wirtschaftliche Zeitreihe in "systematische" Komponenten und 
in eine "Restkomponente" zerlegen läßt. Die inhaltliche Deu- 
tung und relative Abgrenzung der systematischen Komponente 
erfolgt dabei durch substanzwissenschaftlich fundierte Über- 
legungen, während die Restkomponente als diejenige Komponen- 
te verstanden wird, die sich solchen Deutungen entzieht. Ein 
solches Vorgehen läuft zunächst auf das hinaus, was (WALD, 

38) als "äußere" Definition von Komponenten bezeichnet hat. 

Da aber ein derartiges Vorgehen letztlich auf unüberwindli- 
che Schwierigkeiten stößt, wie A. WALD gezeigt hat, und des- 
halb zu keinen realisierbaren Bereinigungsverfahren führt, 
bleibt nur der Ausweg der "inneren" Definition von Komponen- 
ten (WALD, 38) . Dies läuft darauf hinaus, daß man eine ge- 
wisse Anzahl von Eigenschaften der Reihenkomponenten postu- 
liert, von denen man annimmt, daß sie ökonomisch sinnvoll 
oder plausibel sind. Diese Postulate ermöglichen es dann, den 
einzelnen Komponenten je eine gewisse Klasse von Funktionen 
zuzuordnen (z. B. die Klasse der Polynome "niedriger" Ordnung) , 
von denen man der Auffassung ist, daß diese Funktionenklas- 
sen in der Lage sind, eine hinreichend adäquate formale Cha- 
rakterisierung eben dieser aus dem ökonomischen Sachwissen 
heraus formulierten "Eigenschaften" zu liefern. Aber selbst 
wenn alle Konstrukteure von Saisonbereinigungsverfahren 
einheitlich diesen Katalog von Eigenschaften sowie die zu- 
gehörigen Funktionenklassen akzeptierten , so wären trotzdem 
auf dieser Grundlage die verschiedensten Verfahren denkbar, 
da eben diese Grundlage nicht mehr leisten kann, als gewis- 
se Funktionenklassen für jede als ökonomisch relevant er- 
achtete Komponente festzulegen, innerhalb derer dann der 
Statistiker die Möglichkeit hat, eine ihm adäquat erscheinen- 
de Funktion auszuwählen, um damit endgültig die jeweilige 
Komponente zu definieren. Daraus folgt nun: werden z. B. bei 
zwei Verfahren für eine Komponente verschiedene Funktionen aus 
derselben Funktionsklasse ausgewählt, dann ist zu erwarten. 
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daß diese Verfahren zu numerisch u. U. erheblich differieren- 
den Resultaten führen. Diese Unterschiede sind aber ausschließ- 
lich auf verschiedene Definitionen zurückzuführen. Beim Ber- 
liner Verfahren wird z. B. die sog. "glatte" Komponente durch 
ein Polynom dritten Grades definiert, während beim ASA II- 
Verfahren diese Komponente durch eine Funktion definiert wird, 
die sich additiv zusammensetzt aus einem Polynom zweiten Gra- 
des und einer trigonometrischen Funktion mit einer Periode von 
36 Monaten. Selbst wenn sich beide Verfahren nur in der De- 
finition dieser Komponente unterscheiden würden, wären ver- 
schiedene saisonbereinigte Zeitreihen in der Regel die Fol- 
ge. Es wäre aber nicht sehr sinnvoll, darüber zu diskutieren, 
welches Verfahren denn nun eigentlich das "bessere" ist, da 
die unterschiedlichen Resultate eben nur auf eine unterschied- 
liche Komponentendefinition zurückzuführen sind. Eine solche 
Diskussion wäre gleichbedeutend damit, daß man herausfinden 
wollte, welche der beiden Definitionen die "bessere" ist. Es 
ist klar, daß ein solches Unterfangen sinnlos wäre, da die 
Frage, ob eine gewählte Komponentendefinition "besser" oder 
"schlechter" mit der ihr entsprechenden äußeren Komponenten 
übereinstimmt, nicht positiv entschieden werden kann, denn 
dazu müßte ja diese äußere Komponente quantitativ bekannt 
sein, was natürlich nicht der Fall ist, sonst könnte man sich 
ja die ganze Zeitreihenzerlegung sparen (WALD, 38). Da jedes 
Saisonbereinigungsverfahren notwendigerweise eine innere De- 
finition der relevanten Komponenten impliziert, insbesondere 
auch der Saisonkomponenten, sind mit verschiedenen Verfahren 
bereinigte Zeitreihen prinzipiell unvergleichbar, da diese 
Verfahren in der Regel auf unterschiedlichen Definitionen 
beruhen. Der häufig geübte Vergleich alternativ bereinigter 
Reihen, insbesondere von Einzelheiten in ihrem zeitlichen 
Verlauf, entbehrt daher einer sinnvollen Grundlage, und die 
damit verbundenen, oft mit Leidenschaft geführten Kontrover- 
sen, gelten in Wirklichkeit einem Scheinproblem. Uneinge- 
schränkt gilt dies zumindest dann, wenn nicht schlüssig aus 
ökonomischen Überlegungen heraus dargelegt werden kann, daß 
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ein Verfahren tatsächlich "unsinnige" Ergebnisse bringt. Ein 
solcher Nachweis wird jedoch nur ausnahmsweise möglich sein. 

In der Regel basieren derartige Vergleiche auf subjektiven a 
priori-Vor Stellungen darüber, welche Veränderungen einer kon- 
kreten Reihe als saisonal bedingt und welche als nicht-saiso- 

1 ) 

nal bedingt angesehen werden sollen. 

Dieser an sich triviale Sachverhalt kann auch noch von einem 
anderen Aspekt aus betrachtet werden. Es ist heute üblich, 
auch die Analyse ökonomischer Zeitreihen im Rahmen der Theo- 
rie der stochastischen Prozesse zu betrachten. Daß dies mit 
guten Gründen geschieht, ist wohl unstrittig. Eine solche 
Betrachtungsweise impliziert die Vorstellung, daß eine kon- 
krete Zeitreihe eine mögliche Realisation eines hypotheti- 
schen Zufallsprozesses ist, der sich in Teilprozessen zer- 
legen läßt. Der oben skizzierte Ansatz ist nun gleichbedeu- 
tend damit, daß für die Realisationen der Teilprozesse bestimm- 
te Traj ektorieneigenschaf ten postuliert werden. Das läuft 
darauf hinaus, daß unterstellt wird, es sei möglich, diese 
Realisationen durch geschlossene mathematische Ausdrücke 
darzustellen. Offensichtlich läßt sich eine solche Darstel- 
lung auch unter Einbeziehung allen jeweils verfügbaren öko- 
nomischen Sachwissens nicht eindeutig vollziehen. Im übri- 
gen sei hier darauf hingewiesen, daß im Rahmen einer solchen 
Betrachtungsweise offensichtlich nur eine sehr spezielle Klas- 
se von stochastischen Prozessen zugelassen v/ird, nämlich die, 
deren Trajektorien so "brav" verlaufen, daß die erwähnte Dar- 
stellung möglich erscheint. Ob eine derartige Einschränkung 



1) Etwas anders liegen die Verhältnisse, wenn solche Ver- 
gleiche auf der Basis von simulierten Reihen vorgenommen 
werden. Allerdings ist auch hier Vorsicht insofern gebo- 
ten, als darauf zu achten ist, daß diese Reihen nicht so 
konstruiert werden, daß sie a priori ein bestimmtes Ver- 
fahren präferieren, was z. B. im Extremfall dann eindeu- 
tig geschieht, wenn die Komponenten einer künstlichen 
Reihe genau gleich definiert werden wie in einem bestimm- 
ten Säisonbereinigungsverf ahren . 
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auf "Lehrbuch-Prozesse" allerdings der Komplexität ökonomischer 
Prozesse angemessen ist, muß bezweifelt werden. Gibt man diese 
restriktiven Vorstellungen auf, dann sind allerdings auch an- 
dere methodische Konzepte der Zeitreihenanalyse heranzuziehen, 
als die in den heute gängigen Verfahren üblich sind. Darauf 
wird noch ausführlich einzugehen sein. Beiläufig sei hier noch 
darauf hingewiesen, daß die Vorstellung eines "glatten" Ver- 
laufs von Komponenten ökonomischer Zeitreihen möglicherweise 
ihre Wurzel im übermäßigen Gebrauch einfacher gleitender 
Durchschnitte hat, wie er in der traditionellen Zeitreihenana- 
lyse üblich war. Es läßt sich nämlich unschwer nachweisen, daß 
solche Durchschnitte im allgemeinen "schlechte" Tiefpaßfilter 
und eigentlich nur für Komponenten geeignet sind, deren spek- 
trale Masse ausschließlich in der Frequenz Null konzentriert 
ist. (Vgl. z. B. STIER, 37) 

Zweifellos liegt die eigentliche Problematik einer Zerlegung 
von Zeitreihen in Komponenten darin, daß diese Komponenten 
ökonomisch nur vage definiert sind. Man hat nun im Laufe der 
Zeit versucht, diese offensichtliche Schwäche dadurch zu um- 
gehen, daß man stärker formalisierten Komponenten-Def initionen 
den Vorzug gab. Solche ergeben sich mit Hilfe von Konzepten, 
die der Spektralanalyse entlehnt sind. Danach ist z. B. die 
Saisonkomponente als diejenige Komponente einer Zeitreihe de- 
finiert, deren spektrale Masse in den Frequenzen 

2t\]c 

= 1,2,..., 6) konzentriert ist, genauer: innerhalb "klei- 
ner" Intervalle um diese Frequenzen. Sicherlich ist eine der- 
artige Definition von einem formal-mathematischen Standpunkt 
aus betrachtet erheblich präziser als eine nur "deskriptive" 
Definition. Allerdings ist eine solche Definition immer noch 
unvollständig, da ja die Breite dieser Intervalle Undefiniert 
bleibt. Es läßt sich mit Hilfe der Spektralanalyse exakt nur 
definieren, was nicht unter der Saisonkomponente einer ökono- 
mischen Zeitreihe verstanden werden soll. Zwei Extremfälle 
sind denkbar: die Saisonkomponente wird durch ein "Linien- 
Spektrum" definiert bzw. durch ein "flaches Spektrum", d. h. 
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ihre Spektraldichte ist konstant im gesamten Frequenzintervall. 
Eine Definition durch ein Linienspektrum impliziert die Vor- 
stellung, daß saisonale Bewegungen exakt periodisch verlaufen, 
eine Vorstellung, der keinerlei reale Bedeutung zukommt, da 
exakt-periodische Bewegungsabläufe mathematische Fiktionen dar- 
stellen, die sich selbst im Bereich der Naturwissenschaften, 
wie z. B. in der Astronomie, als unrealistisch erweisen. An- 
dererseits würde eine Definition durch ein flaches Spektrum 
bedeuten, daß saisonale Bewegungen völlig unstrukturiert ver- 
laufen, also einen reinen Zufallsprozeß darstellen. Diese Vor- 
stellung ist natürlich absurd, denn sie wäre gleichbedeutend 
damit, daß das Phänomen "Saison" realiter gar nicht existier- 
te. Realistisch ist deshalb nur der Fall, der zwischen die- 
sen beiden Extremen liegt. Je "regelmäßiger" eine Saisonkom- 
ponente verlaufen soll, um so enger sind die Intervallbreiten 

zu definieren und umgekehrt. Aber so lange diese Breiten nicht 

1 ) 

objektiv festzulegen sind , bleibt eben eine auf der Grund- 
lage der Spektralanalyse vollzogene Definition der Saisonkom- 
ponente, d. h. also, eine Definition im Frequenzbereich, un- 
vollständig . 

Es ist heute üblich geworden, die Transferfunktionen der bei 
Saisonbereinigungsverfahren eingesetzten Filter zu berechnen 
und explizit anzugeben. Diese Transferfunktionen lassen letz- 
ten Endes die im jeweiligen Verfahren implizierte Definition 

2 ) 

der relevanten Komponenten im Frequenzbereich erkennen , 
sofern sie überhaupt interpretierbar sind. Das bedeutet aber 



1) Eine Untersuchung dieser Möglichkeit erfolgt in Abschnitt 

3.2.3. 

2) Den heute praktizierten Verfahren ist gemeinsam, daß zu- 
erst Filter konzipiert und dann deren Transferfunktionen 
analysiert werden. Zweckmäßiger scheint der umgekehrte 
Weg zu sein. Darauf wird noch in Abschnitt 3.3 eingegan- 
gen . 
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1 ) 

wiederum, daß z. B. unterschiedliche Bandbreiten bei der 
saisonalen Komponente eine unterschiedliche Definition die- 
ser Komponente implizieren. Somit ergibt sich also auch bei 
einer Betrachtung der Definitionsprobleme im Frequenzbereich, 
daß die mit Hilfe von verschiedenen Verfahren bereinigten 
Zeitreihen prinzipiell unvergleichbar sind. 



1 ) Auf eine Diskussion des in der Literatur mit verschiedenen 
Inhalten belegten Begriffes "Bandbreite" sei hier verzich- 
tet. Hier ist einfach die Länge der oben erwähnten Inter- 
valle gemeint. Vorausgesetzt ist dabei nur, daß sich die- 
se bei den praktizierten Verfahren hinreichend genau er- 
mitteln und vergleichen lassen. 
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3.2.3. Verknüpfung von Zeitbereichsanalyse und Frecruenz- 
bereichsanalyse? 

Bei den Saisonbereinigungsverfahren, die von der Annahme aus- 
gehen, daß sich die Trajektorien der einzelnen Teilprozesse 
durch geschlossene mathematische Ausdrücke mindestens inner- 
halb im Vergleich zur Gesamtlänge der vorliegenden Zeitreihe 
relativ kurzen Zeitintervalle darstellen lassen, stellt sich 
das Problem, die unbekannten Parameter, die in diesen Aus- 
drücken enthalten sind, zu schätzen. Sowohl beim Berliner 
Verfahren als auch beim ASA-II-Verf ahren z. B. werden dafür 
Schätzfunktionen in Form von gleitenden Durchschnitten ver- 
wendet, deren Gewichte mit Hilfe der Methode der kleinsten 
Quadrate bestimmt werden. Dieser Ansatz, der in seiner ur- 
sprünglichen Form auf WEICHSELBERGER zurückgeht, führt zu 
linearen erwartungstreuen Schätzfunktionen mit minimaler Va- 
rianz (WEICHSELBERGER, 40) , Somit genügen diese gleitenden 
Durchschnitte zunächst einmal Optimalitätskriterien , deren 
Beachtung im Rahmen der statistischen Schätztheorie üblich 
ist. Allerdings sind diese Kriterien nur im Zeitbereich de- 
finiert. Deshalb sind so konstruierte gleitende Durchschnit- 
te oder Filter nur in diesem Bereich als "optimal" anzuse- 
hen. Damit ist aber noch nichts über irgendwelche "Optimal- 
Eigenschaf ten" der Transfer Funktionen dieser Filter im Fre- 
quenzbereich ausgesagt. Auffallend ist bei den erwähnten 
Verfahren die Tatsache, daß solche weder explizit formuliert 
noch implizit in den verwendeten Algorithmen erkennbar sind. 
Vielmehr werden die Transf erfunktionen der jeweiligen Filter 
nur nach Augenschein auf unerwünschte Absorptions- und Ver- 
stärkungseffekte untersucht und schließlich diejenigen 
Filter herangezogen, deren Transferfunktionen jeweils "am 
besten" mit den intendierten "idealen" Transferfunktionen 
überein st immen . Es dürfte klar sein, daß ein solches Verfah- 
ren unbefriedigend ist, da bei der Filterauswahl weder Ap- 
proximations- noch sonstige Gütekriterien exakt definiert 
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und berücksichtigt werden. Infolgedessen ist die Filterauswahl 
in diesen Verfahren letzten Endes das Resultat eines mehr oder 
weniger planlosen Herumprobierens. Nur so ist es zu verstehen, 
daß z. B. beim Berliner Verfahren über 200 Filter auf ihre 
Transferfunktionen hin untersucht wurden (NULLAU, 24) . An- 
scheinend ist diese "Probiererei" bei diesem Verfahren bis 
heute nicht zu Ende gekommen (Vgl. NOURNEY, 23) . Bei einem 
derartigen Vorgehen steht im Hintergrund die Vorstellung 
einer "idealen" Transferfunktion. Diese Vorstellung impli- 
ziert aber bzw. ist äquivalent einer Definition der jeweili- 
gen Komponente im Frequenzbereich. Wenn z. B. als ideale 
Transferfunktion eines Saisonkomponenten-Filters die Funktion 



T (A) 



r 1 für X = k X + 6 
0 sonst 



A 



2tt 

12 ' 



k = 1 ,2, ... ,6 



angesehen wird (Vgl. z. B. HEILER, 13), dann heißt das doch 
nichts anderes, als daß damit eine Definition dieser Kompo- 
nente im Frequenzbereich formuliert wird. Konkret ist in die- 
sem Beispiel die Saisonkomponente als Realisation eines sto- 
chastischen Prozesses definiert, dessen spektrale Masse in 
den angegebenen Frequenzbändern der Länge 2 6 konzentriert 
ist. Wie oben schon ausgeführt, ist eine so vollzogene Defi- 
nition aber so lange unvollständig, so lange die Größe 6 
nicht bestimmt ist. Die Frage ist nun hier, ob sich diese 
Größe prinzipiell für diesen und/oder andere Idealfilter in 
irgendeiner sinnvollen Weise, d. h. also mittels Kriterien 
der statistischen Schätztheorie, objektiv festlegen läßt. Wä- 
re das möglich, dann ließe sich in der Tat eine Brücke schla- 
gen zwischen einer Zeitreihenanalyse im Zeitbereich und einer 
Zeitreihenanalyse im Frequenzbereich, also Bereichen, die bis- 
her offensichtlich nur lose miteinander verknüpft sind. 

Die Untersuchung dieses Problems erfordert recht umfangreiche 
Ableitungen. Dabei werden eine Reihe von idealen Vorausset- 
zungen gemacht, was für eine Analyse eines prinzipiellen Pro- 
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blems sicherlich gerechtfertigt ist. Falls eine positive Pro- 
blemlösung existiert, müßte dann in einem zv/eiten Schritt nach 
der praktischen Bedeutung bzw. Verwertbarkeit dieser Lösung 
gefragt werden . 

Ausgangspunkt ist zunächst ein denkbar einfaches Zeitreihen- 
modell: 



X(t) = S(t) + N(t) 



Dabei sei S(t) ein Saisonprozeß und N (t) ein Restprozeß. Es 
wird nun angenommen, daß die analytische Form dieser Prozesse 
bzw. ihre Spektren bekannt sind. Für S(t) wird der allgemein 
bekannte Hannansche Saisonprozeß unterstellt und für N(t) ein 
"white-noise"-Prozeß . Für diesen Saisonprozeß ergeben sich 
für Autokovarianzfunktion, Varianz und Spektrum die Ausdrücke: 



Y (t) 






6 

Var S(t)=^|^ 




COST Y . 

' 1 



» |Pil<i 



T — 0, +1 +2,... 



und 



S* (X) 

SS 



— - E 
4tt i=1 




o 



2 

i 



2 cos(X-X^) 



2 
a . 

1 

1+p^ - 2 p^ cos(X+ X^) 



mit 



2tt± 

12 



i=1,2,...,6 , |X| <TT 



Für die Koeffizienten (t) und 3j^(t) des trigonometrischen 
Polynoms gelten dabei die Beziehungen 
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cti(t) = p^a^(t-1) + £^(t) 

B^(t) = p^ß^(t-1) + ^^(t) 

E[a^(t) ]=E[3^(t) ] = 0 , 

E[e^ (t) ]=E[C^ (t) ] = 0 

wobei p. der Autokorrelationskoeffizient erster Ordnung von 
1 2 

a^(t) bzw. B^(t) ist. bezeichnet die Varianz von e^(t) bzw 
C^(t). (Vgl. HANNAN, 12). 

Da N(t) als reiner Zufallsprozeß mit Erwartungswert Null ange 
sehen wird, gilt für das Spektrum dieses Prozesses 

S*n(^) = 2 ^ N(t) = const. , |X| < TT 

Da die Prozesse S (t) und N (t) schwach-stationär sind, ist 
X (t) schwach-stationär. Für den Erwartungswert, die Varianz 
und das Spektrum dieses Prozesses gelten die Ausdrücke: 

E[X(t)] = E[S(t)] + E[N(t)] = 0 
Var X(t)= Var S(t)+ Var N(t) 



S*^(A) = 









|A I < TT 



wenn zusätzlich noch postuliert wird: 



Cov[S (t) ,N (t) ] = 0 , V t 



d. h. wenn also die Prozesse S(t) und N(t) als unkorreliert 
angenommen werden. 

Für das Folgende ist es zweckmäßig, mit einem normierten Spek 
trum für den Prozeß X (t) zu arbeiten. Dieses ergibt sich aus 
S* (X) durch Division mit Var X(t), wobei gilt 

= } S*^(X) d A 

-n 



Var X (t) 
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Für die Spektren der Teilprozesse S(t) und N (t) ergeben sich 
dann die Ausdrücke 





- S|g(A) 


/ 


Var 


X(t) 




S*^(A) 


/ 


Var 


X(t) 



und es gilt für das normierte Spektrum des Gesamtprozesses 



Sxx<^) = Sgg(A) + S^^(A) , |A|< TT 



Nun ist 



Var N(t) _ Var S(T) 

Var X (t) Var X (t) ^ 



Daraus folgt 



1 - 2^ S^^(A) < 0 



d. h. 



NN 2 TT 



Weiter folgt aus 



N 



0 • 



^NN 



1 Var N(t) 
2tt Var X (t) 



die Beziehung 



N ^ _1 Var N(t) 

0 2tt Var S (t) + Var N (t) 

_± 1 

2tt Var S(t) 

Var N(t) 



Für das vorliegende Zeitreihenmodell, das offensichtlich vom 
Typ "Signal plus noise" ist, muß sinnvollerweise angenommen 
werden, daß gilt 
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Var S(t) >> var N(t) 

Somit gilt für das Spektrum des Rauschprozesses 

O < Nq << 1 

Ziel der Zeitreihenanalyse sei es nun, das "Signal" S(t) durch 
einen linearen Filter zu schätzen. Die Schätzfunktion sei 

00 

S(t) := Z hg(k) X(t-k) 

k=-oo 

mit E[S (t) ] = 0 

wobei hg (k) die Impuls-Antwortfunktion des Filters ist. (Das 
Suffix S soll auf "Saison" hinweisen) . Vorausgesetzt wird, 
daß die Zeitreihe X (t) "unendlich" lang ist. Der Schätzfehler 
sei durch die Differenz 

F(t) := S(t) - S(t) 



definiert . 

Der Filter soll nun so gewählt werden, daß gilt 

Q := E[F^(t)] = E[S(t) - S(t)]^ = Min. 

d. h. S (t) soll eine lineare Schätzfunktion für S (t) sein mit 
minimalem mittlerem quadratischen Fehler. Existiert ein sol- 
cher Filter, dann ist 

N(t) := X(t) - S(t) 

eine Schätzfunktion für N (t) mit minimalem mittlerem quadra- 
tischen Fehler, denn es ist 

E[N(t) - N(t)]^ = E[N(t) - {X(t) - S(t)}]^ 

= E[N(t) - {S(t) + N(t) - S(t)}]^ 

- S(t) 



= E[S(t) 
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Umgekehrt folgt daraus, daß 

S(t) = X(t) - N(t) 

eine Schätzfunktion für S (t) ist mit minimalem mittlerem qua-' 
dratischen Fehler, wenn 

^ oo 

N(t) := Z h^(k) x(t-k) 

k=-oo 

eine solche ist für N(t), wobei h^ (k) die Impuls-Antwortfunk- 
tion des Noise-Filters ist. In diesem einfachen additiven Zeit- 
reihenmodell ist es also prinzipiell gleichgültig, für welchen 
Teilprozeß ein Filter konzipiert wird, die jeweils interessie- 
rende Komponente ergibt sich immer durch eine einfache Ope- 
ration. Der im Folgenden öfters verwendete Begriff "Saison- 
filter" ist doppeldeutig. "Filtern" kann hier bedeuten: Schät- 
zen von S(t) oder Eliminierung (d. h. also "Ausfiltern") von 
S(t), also primär Schätzen von N (t) und erst sekundär Schät- 
zen von S (t) durch Differenzenbildung. 

Aus der Autokovarianzfunktion 

Rpp (t) := E[S(t) -S (t) ] [S (t+T) - S (t+i) ] , t = 0,+1 ,+2, . . . 

des "Fehler-Prozesses" F (t) ergibt sich für den mittleren 
quadratischen Fehler 



Q = Rpp(O) 



Nun ist aber 



1 

= 2? Sff(M d X , I X I < TT 

-TT 

wobei Spp(A) das Spektrum des Fehlerprozesses F(t) ist. Kann 
Spp(A) bestimmt werden, dann ist der mittlere quadratische 
Fehler des Saisonfilters S (t) im Freqenzbereich gegeben und 
die oben diskutierte Problemstellung kann untersucht werden. 
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Die Autokovarianzfunktion R„„(t) des Fehlerprozesses F (t) läßt 

r r 

sich nun als Summe von Autokovarianz- bzw. Kreuzkovarianzfunk- 
tionen dar stellen: 



Rpp(T) = Rgg(x) - Rgg(i) - Rgg(T) + Rgg(T) > t=0,+1,+2,. 



Rgg(T) := E[S(t) S(t+T)] , T = 0, +1, +2,... 

Rgg(T) := E[S(t) S(t+x)] , X = 0, +1, +2,... 

Rgg(x) := E[S(t) S(t+x)] , X = 0, +1, +2,... 

Rgg(x) := E[S(t) S(t+x)] , X = 0, +1, +2,... 

“i At 

Durch Multiplikation von Rpp ( t) mit e und Summation er- 

gibt sich: 



SppCA) = Sgg(X) - Sg^CA) - Sgg(A) + Sgg(A) , |A| < ^ 

(Da sowohl S (t) als auch S(t) schwach stationär sind, existie- 
ren Sgg(A), Sgg(A), Sgg(A) sowie Sgg(A) ). 

Bekannt ist zunächst nur Sgg(A) . Bestimmt werden müssen noch 
das Spektrum Sgg(A) sowie die Kreuzspektren Sgg(A) und Sgg(A). 

Es ist 

Sgg(A) := Z e , x = 0,+1,+2,... 



= Z E[S (t) S (t+x) ]e 

q; = — oo 



S(t+x) = E h^(k) X(t+T-k) , T = 0,+1,+2,... 

k=-oo ^ 
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Das Signal (t) läßt sich ebenfalls formal durch eine "Filter- 
gleichung" darstellen: 

oo 

S(t) = Z h(k) S(t-k) 

k=-oo 

wenn für die Impul s-Antwortfunktion h(k) gilt: 



h(k) 



r1 für k = O 
für k + 0 



Dies ist nämlich die impuls”Antwortfunktion eines "Allpass- 
Filters, d. h. also eines Filters mit der Frequenz -Antwort- 
funktion 



T(A) = 1 , |X |<TT 

denn es ist: 

1 ^ i X k 

h(k) = 2^ f T(X) e^^^ dX 



und daraus folgen 



sowie 



d. h. 



h(k) = 1 , k = 0 
1 

h(k) = sinTTk, 
h(k) = 0 , 



k 4 0 
k 4 O 



Für das Kreuzspektrum zwischen S(t) und S(t) folgt damit der 
Ausdruck: 



Sgg (A) 



oo 

z 

T=* 



oo oo 

E[ Z h(k)S(t-k) Z 
k=-od k=- 



hg (k) X (ti-T-k) ]e 



oo oo oo 

= Z h(k) Z hg(l) Z Rg^(T+k-l) e 

k=— oo 2^=— oo — oo 



(T+k-1) := E[S(t-k)X(t+T-l) ] 



mit 
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Mit 



m := T + k - 1 



ergibt sich daraus 






^ ^ ^ ^ ^ r, , N (m+l-k) 

E h(k) E hg(l) E Rg^ (m) e 



°° -i Ak “ 

= E h(k) ^ E 

k=— 1= 



hs(l)e 



-iAl 



Rsx (in) 



-iXm 

e 



= T(X) Tg(X) Sg^(X) 

Dabei sind T(X) die Kon jugi er t- Komplexe der Frequenz -Antwort - 
funktion des Allpass-Filters, Tg(X) die Frequenz-Antwortfunk- 
tion des Filters S(t) und Sg^(X) das Kreuzspektrum der Prozes- 
se X(t) und S(t). Da T(X) reell ist und für alle Frequenzen den 
Wert Eins annimmt, kann schließlich geschrieben werden: 

^SS ^ ^SX 

Ganz analog erhält man ähnliche Beziehungen für Sgg(X) und 
Sgg(X). Auf eine explizite Herleitung kann deshalb verzich- 
tet werden . Es ist: 

^SS " '^S ^XS 

und 

= l^sl" ' 

wobei die Kon jugiert-Komplexe der Frequenz-Antwortfunktion 
des Saisonfilters S(t), |Tg| die Transferfunktion dieses Fil- 
ters und S „{X) das Kreuzspektrum der Prozesse X(t) und S (t) 

X o 

bedeuten . 
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Mit diesen Ausdrücken läßt sich dann für das Spektrum S„„(X) 

r r 

des Fehlerprozesses F (t) schreiben: 

Spp(X) = Sgg(X) - Tg(A)Sgjj(X) ~ T g ( A ) S g ( A ) + | T g ( A ) | ^ S ( A ) 

Dieser Ausdruck läßt sich noch vereinfachen. Es ist nämlich: 

Rg^(x) = E[S(t)S(t+x) ] 

= E{S(t) [S(t+x)+N(t+x) ] 

= E[S(t)S(t+x) ]+E[S(t)N(t+x) ] 

= «SS ' 

da S(t) und N (t) voraussetzungsgemäß unkorreliert sind. Da- 
her ist 

Ebenso einfach folgt die Beziehung 

= ^ss(^) 

Zusammen mit dem Ausdruck 
folgt dann 

l^+ t1-2 Re Tg(A) + |Tg(A) |^] 

wobei Re Tg(A) den Realteil der Frequenz-Antv.^ortf unktion des 
Saisonfilters bezeichnet. 

Damit ergibt sich schließlich für den mittleren quadratischen 
Fehler der Schätzfunktion S(t), d. h. also des Saisonfilters, 
der Ausdruck: 

2 I ^ '^S (^)<3A+^ / Sgg(A) [Tg(A)-1]2 dA 

— TT -TT 
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wenn Tg(A) als rein reell angenommen wird, was für die nach- 
stehenden Untersuchungen zweckmäßig ist, da in vorliegendem 
Zusammenhang evtl. Phasenverschiebungen des Saisonfilters 
nicht interessieren. Mit diesem Ausdruck ist eine Beziehung 
zwischen dem Zeitbereich und dem Frequenzbereich hergestellt. 

Es zeigt direkt die Abhängigkeit des mittleren quadratischen 
Fehlers des im Zeitbereich definierten Filters S(t) von sei- 
ner Frequenz-Antwortfunktion bzw. Transferfunktion. 

Die nachfolgenden Untersuchungen basieren nun auf folgenden 
Überlegungen: Ausgangspunkt ist die Wahl verschiedener "idea- 
ler" Transf erf unktionen , die jeweils bis auf einen zu bestim- 
menden Parameter spezifiziert sind. Wie oben dargelegt wurde, 
ist mit der vollständigen Spezifikation einer Transferfunk- 
tion die jeweils auszufilternde Komponente, also hier: die 
Saisonkomponente, im Frequenzbereich eindeutig definiert. 

Es wird nun versucht, den a priori unbestimmten Parameter 
jeweils "optimal" festzulegen. Dabei soll immer derjenige 
Parameterwert als optimal angesehen werden, der den mittle- 
ren quadratischen Fehler der Schätzfunktion S(t) minimiert. 
Existieren derartige optimale Transferfunktionen, dann könn- 
te prinzipiell die oben diskutierte Unvollständigkeit von 
Komponentendefinitionen im Frequenzbereich behoben werden. 

Für den Problemkreis "Saisonbereinigung" hieße das, daß die 
relevanten Bandbreiten eindeutig und objektiv bestimmbar wä- 
ren. Damit könnte die oben als bislang unvollständig charak- 
terisierte Definition der Saisonkomponente einer ökonomischen 
Zeitreihe durch Konzepte aus der Spektralanalyse nun ein- 
deutig vollzogen werden. 

Im einzelnen werden im Folgenden fünf Transf erfunktionen unter- 
sucht, die immer wieder in der Literatur und/oder bei Diskus- 
sionen als "ideal" für die Saisonbereinigung angesehen werden. 
Außerdem liegt die Wahl dieser Transf er funktionen hier nahe, 
da bei der im Zeitreihenmodell zugrunde gelegten Hannanschen 
Saisonfigur das korrespondierende Spektrum S*g(A) die Eigen- 
schaft hat, daß seine spektrale Masse in den saisonalen Fre- 
quenzen konzentriert ist. 
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Von den zu analysierenden fünf Transferfunktionen sind drei 
Filtern zuzuordnen, die als "Selektionsfilter" bekannt sind. 
Die Transferfunktionen solcher Filter haben die Eigenschaft, 
daß sie intervallweise konstant sind. Diese sollen zuerst be 
trachtet werden. Anschließend werden zwei Transferfunktionen 
eines anderen Filtertyps untersucht. Prinzipiell lassen sich 
die den einzelnen Transferfunktionen korrespondierenden 
Filterkoeffizienten durch die Beziehung 

1 ^ i A k 

hs(k) = 2¥ / Tg(A) e'^^'^dA, k = 0,+1 , +2,.. 



ermitteln. Darauf wird aber im Folgenden verzichtet, da es 
im vorliegenden Zusammenhang nicht um das Problem der Kon- 
struktion von Filtern geht. Dazu soll erst weiter unten Stel 
lung genommen werden. Es sei hier lediglich schon angemerkt, 
daß der mittlere quadratische Fehler von S(t) auch dann von 
Null verschieden wäre, wenn die zu diskutierenden Ideal- 
Filter realisiert werden könnten. Ideal-Filter implizieren 
also nicht eine fehlerlose Schätzung, wie das Adjektiv 
"ideal" eventuell suggerieren könnte. Das wäre im vorliegen- 
den Zeitreihenmodell nämlich nur dann der Fall, wenn sich 
die Spektren von S(t) und N(t) nicht überlappen würden, was 
offensichtlich der Fall ist. 

Die erste Transferfunktion ist wie folgt definiert: 



f O für 0 < A < 



2tt 

12 



Tg,I(A) 






1 f ür - 6 < A < ^ + 6 , k = 1 , 



12 



O für ^ + 6 < A< ^2 



12 

2tt (k+1 ) 



<5 , k = 1 , 



' 1 für 



6 < A < 



,5 

,5 
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12 12 12 12 12 



Dies ist die Transfer funktion eines multiplen Bandpassfilters 
mit der Bandbreite 2 6. Diese Funktion ist vollständig spezi- 
fiziert, wenn der Parameter 6 konkretisiert ist. Damit wäre 
dann die Saisonkomponente im Frequenzbereich vollständig de- 
finiert. Die Frage ist nun, ob sich 6 im oben diskutierten 
Sinn bestimmen läßt. 

Da für die folgenden Abteilungen immer gilt: 



^NN * ^ ^ ^NN ^ ^ ^ 






Tg (-A) = Tg (X) 



I X I < TT 
|X I < TT 

I X I < n 



kann für den mittleren quadratischen Fehler in Abhängigkeit 
von der jeweils gewählten Transferfunktion geschrieben wer- 
den : 
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TT 

Tg(X) dA + l/Sgg(A) [Tg(A)-1 ]^dA 



0 < X < TT 



Wie aus der Definition von ersichtlich ist, muß für 

0,1 



die Berechnung von Q [T„ ^ {^) ] 

0,1 



das Intervall [0,tt] in 12 
Teilintervalle zerlegt werden. Im einzelnen ergeben sich für 
diese Intervalle folgende Resultate: 



a) Intervall J 



O 



:=[ 0 ,|| - 6 ) 



2tt 

l12 
^ " 0 



- 6 

Sgs(A) dA.= 



1 



TT VarX(t) 



2tt . 

T2 " ^ 

/ S* (A) dA 
O 



Grundsätzliche Schwierigkeiten bereiten bei dieser wie auch 
bei den anderen Transferfunktionen die jeweiligen Integrale 
der Funktion S*g(A). Diese führen nämlich zu ziemlich kompli- 
zierten Ausdrücken, so daß die Frage, ob für Q[T„(A)] ein 

O 

eindeutiges Minimum bezüglich des Parameters 6 existiert, 
nicht mehr auf analytischem Weg beantwortet werden kann . 

Ein Ausweg aus dieser Problematik kann wie folgt gefunden 
werden: da S*g(A) stetig ist im Intervall [0,tt] und damit 
auch in jedem Teilintervall von [0,tt] gilt nach dem Mittel- 
wertsatz der Integralrechnung für obiges bestimmtes Inte- 
gral: 



2tt 

12 



- 6 



,2tt 



/ S*3(X) ax = (^ - 5) S*3(5„) 



0 < < 



2tt 

12 



- 6 



Deshalb ist; 



^ ^ss 

TT Var S (t) 



f 2n 

M2 



6 ) 



Jq = [0,fl - 6) 
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Bei den anderen Intervallen wird analog vorgegangen: 



b) Intervalle J,_ := - S, +6), k = 



’k 

2irk 
1 ^2 

1 = - / 
•^k 2 tiJc 

12 



+ 6 



Nq dX 



2 

- 6 
TT 0 



- 6 



c) Intervalle J* :=(^ + 6, k=1,...,5 



2tt (k+1) 

1 2 ^ 

" TT VarX(t) 



2it (k+1 ) 



2jrk 

12 



t 6 



12 

/ 

2TTk 

12 



Ss*(A)dA 



+ 6 



1 ,2. 



TT VarX (t) M 2 



f^n 2uk , X ^ r / 2TT(k+1) . 

(tT ” 26) , ^ ^ ^k f2 ^ ^ 



d) Intervall •= [tt-6 , tt ] 



^^2 tt-6 ° 



- 6 

TT 0 



Damit ergibt sich 



Q[T (X),6] = Q + 1 + £ Q,*+Q-r 

'^0 k=1 ^k k=1 -^k '^12 



d. h. 

Q[Ts^i(X),63= 1[11Nq- fS3*(Co)+2j S3*(C^)}] + 

I 



^ 6 V^rX(t) + ,^"ss(5k)] 

k=1 



^ - 6 



hLJS ^ ^ . r . 2 TT (k+ 1) 
12 ^ ^k 



12 



- 6 
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Der mittlere quadratische Fehler des multiplen Bandpassfil- 
ters ist also eine lineare Funktion seiner halben Bandbrei- 
te 6. Deshalb existiert kein Minimum für Q[T„ ^(X)] bezüg- 
lieh dieses Parameters. 

Die zweite Transf erf unktion hat folgende Gestalt: 



(^ für O < A < - 5 






1^ 27Tk r ^ ^ 27Tk . . 

|0 für — - 6 < X < —^2 + <5 



!1 für ^ + 6 < 1 < 



k= 1 , . . . , 5 
- 6 , k=l , . . . ,5 



.0 für TT - 6 < X < TT 




12 12 12 12 12 

Diese multiple Bandsperre eliminiert, wie oben kommentiert, 
die Saisonkomponente und stellt deshalb zunächst ein Filter 
zur Ermittlung von N(t) dar. 

Die Berechnung von Q[T„ (X)] verläuft ganz analog zur Be- 

o ^ X -L 

rechnung von Q[T^ (X)]. Deshalb sei auf eine explizite Ab- 

0,1 

leitung verzichtet. Dasselbe gelte auch für die restlichen 
Transferfunktionen . 
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Es ist: 






{2 



5 

E 

k=1 






27Tk 

12 



6 < 



27Tk 

12 



6 



k= 1 , . . . , 5 



Auch für diesen Filter existiert kein minimaler mittlerer qua- 
dratischer Fehler bezüglich des Parameters 6 . 

Die dritte Transferfunktion stelle eine Modifikation der zwei- 
ten Transferfunktion dar. Sie nimmt in den Sperrbereichen nicht 
den Wert Null an, sondern ist dort jeweils gleich N^, d. h. der 
entsprechende Filter soll möglichst wenig die Noise-Komponen- 
te tangieren. Diese Forderung wird in den einschlägigen Dis- 
kussionen immer wieder erhoben. 

Diese Transferfunktion ist also: 



k=1,...,5 
6 , k=1 , . . . , 5 



1 für O < A < 



12 






27Tk X/A/2TTk,.. 

'^O “12 ■ ^ ^ ^ ^ ^ ^ * 



2TTk , r ^ A / 2it (k+1 ) 
1 für — + 6 < A < Y 2 — 



\^Q fürTT-6<A<7T 
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Transferfunktion III 




^ ^ ^ §JL 10TT TT A 

12 12 12 12 12 



und für den zugehörigen mittleren quadratischen Fehler erhält 
man den Ausdruck 



1 

^*•'^3,111 ^ "VarX (t) 



, ^SS ■*■^33 1 2^ ^ ^ "^0 



27Tk 

12 



6 < ^k < 



2TTk . 

12 



6 



TT - 6 < 5 ^ 2 < 'fl' 






Auch für diesen Filter ergibt sich wiederum eine lineare Be- 
ziehung zwischen Schätzfehler und Bandbreite. 

Bei der vierten Transf erfunktion handelt es sich um die Trans- 
ferfunktion eines multiplen Kerbenfilters. Ein solcher Filter 
wird z. B. von (MEISNER, 22) als "ideal" für die Saisonbe- 
reinigung angesehen. Die Transferfunktion lautet: 




1 



2-jjk , . ^ , 2tf (k+1 ) 

für yy- + 6<A< — ^ 



6 , k=1,. 




Offensichtlich eliminiert dieser Filter den Saisonprozeß 
Für diesen Filter ergibt sich: 
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Auch hier zeigt sich wieder ein linearer Zusammenhang zwi- 
schen dem Parameter 6 und dem mittleren quadratischen Feh- 
ler des Filters. Dasselbe ergibt sich auch für die letzte 
zu betrachtende Transferfunktion, die eine Modifikation der 
Transferfunktion IV dar stellt, wobei diese Modifikation aus 
demselben Motiv wie bei der Transferfunktion II vorgenom- 
men wurde. 

Für diese Transferfunktion gilt: 




2tt 4tt 6tt 8tt 1 Ott tt A 

12 12 12 12 12 
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Es ist: 

Q[Tg^^(A) ,6] 



111 1 
-[-^No(No-1)+ 6 VarX(t) ^ 



{3ii^+9-8^, } X 



0 






2TTk 
1 2 



6 < 



2TTk 

12 



k=1 , . . . ,6 



2TTk 

12 



2TTk 

12 



6 



k=1,...,5 



Somit ergibt sich in allen fünf Fällen, daß der Parameter 6 
nicht optimal festzulegen ist. Dieses sicherlich etwas über- 
raschende Resultat könnte nun darauf zurückzuführen sein, 
daß das zugrundeliegende Zeitreihenmodell zu einfach ist. 

Das Modell kann nun in zweifacher Hinsicht erweitert werden. 
Zunächst soll der Prozeß N(t) nicht mehr ein "white-noise"- 
Prozeßt sein, sondern ein Prozeß, dessen Spektrum nicht 
flach ist. Das Spektrum eines solchen "nicht-saisonalen" 

Prozesses könne in seinem Verlauf in Anlehnung etwa an 

1 ) 

(GRANGER, 10) als monoton fallend angenommen werden. 

Es ist jedoch leicht einzusehen, daß auch für dieses erwei- 
terte Zeitreihenmodell die oben abgeleitete Linearität der 
Fehlerfunktionen Q[Tg(A)] bestehen bleibt. Offensichtlich 
wird von dieser Modifikation nur der erste Integralausdruck 
berührt, da jetzt nicht mehr konstant ist. Mit Hilfe 

des Mittelwert Satzes der Integralrechnung ergibt sich je- 

2 

doch wiederum nur ein Integral über Tg(A) und damit ändert 
sich prinzipiell nichts an den obigen Ergebnissen. Es sei 
noch angemerkt, daß für dieses erweiterte Zeitreihenmodell 



1) Die Frage, ob GRANGER ' S "typical spectral shape of an 
economic variable", das Resultat der verwendeten Schätz- 
technik bzw. spezieller Filteroperationen ist, soll hier 
nicht diskutiert werden, da es auf derartige Subtilitäten 
für die vorliegende Problemstellung nicht an komm t . 
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die Transferfunktionen III und V offenbar nicht mehr sinnvoll 
sind. Die zweite Erweiterung bezieht sich auf die oben getrof- 
fene Voraussetzung der Unkorreliertheit der beiden Teilprozes- 
se S(t) und N(t) . Läßt man diese Voraussetzungen fallen, wobei 
jetzt N(t) wieder der obige nicht-saisonale Teilprozeß mit 
nicht flachem Spektrum ist, dann wird der Ausdruck Q[Tg(X)] 
wesentlich komplizierter, da das Kreuz -Spektrum zwischen bei- 
den Teilprozessen zu berücksichtigen ist. Es erscheint aber 
prinzipiell unmöglich zu sein, einen analytischen Ausdruck 
für dieses Kreuz -Spektrum anzugeben. Während ein im wesentli- 
chen monoton fallender Verlauf des Spektrums von N(t) als 
empirisch gesichert gelten kann, und die hier zugrundegeleg- 
te Hannansche Saisonfigur aufgrund ökonomischer Überlegungen 
mindestens plausibel erscheint, können derartige allgemeine 
Aussagen bezüglich des Kreuz-Spektrums der beiden Teilprozes- 
se nicht gemacht werden. Insofern erscheint eine Ableitung der 
entsprechenden Relationen für dieses erweiterte Modell über- 
flüssig. 

Somit ergibt sich das Resultat, daß die angestrebte Verbin- 
dung zwischen Zeit- und Frequenzbereichsanalyse mit Hilfe 
der oben diskutierten Idealfilter nicht möglich ist. Dies 
impliziert insbesonders , daß mit ihnen eine objektive Ver- 
vollständigung der Definition der Saisonkomponenten einer 
ökonomischen Zeitreihe unter Zugrundelegung von Kriterien 
der statistischen Schätztheorie nicht zu erreichen ist. Die- 
ses Ergebnis wurde unter bestimmten Voraussetzungen abgelei- 
tet: Spezielle Ideal-Filter wurden gewählt und außerdem an- 
genommen, daß der saisonale Prozeß ein spezielles und be- 
kanntes Spektrum besitzt. Sicherlich sind andere Filterty- 
pen und andere saisonale Spektren denkbar. Hier kam es nur 
darauf an zu zeigen, daß das gestellte Problem selbst unter 
diesen relativ einfachen Voraussetzungen nicht lösbar ist. 

Dies läßt es mehr als zweifelhaft erscheinen, ob komplizier- 
tere Voraussetzungen (andere Filter, unbekanntes saisonales 
Spektrum) eine Lösung bringen, d. h. ob es in praxi möglich 
sein wird, eine im Sinne der statistischen Schätztheorie 
objektivierte verfahrensimmanente Definition dessen anzuge- 
ben, was üblicherweise vage mit dem Terminus "Saisonkompo- 
nente" bezeichnet wird. 
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3.3. Saisonbereinigung und rekursive 
Filter 

Die oben erzielten Ergebnisse lassen eine objektive Verknüp- 
fung von Zeitreihenanalysen im Zeitbereich und im Frequenzbe- 
reich als unmöglich erscheinen. Somit existieren für den Pro-- 
blemkreis "Saisonbereinigung" grundsätzlich zwei verschiede- 
ne Zugangsmöglichkeiten. Die erste Möglichkeit besteht darin, 
daß man von der Hypothese ausgeht, eine ökonomische Zeitrei- 
he lasse sich in Komponenten zerlegen, wobei traditioneller- 
weise von drei bzw. vier derartigen Komponenten ausgegangen 
wird. Darüber hinaus sind weitere Hypothesen erforderlich, 
die sich auf die Art der Verknüpfung sowie die Trajektorien 
dieser Komponentenprozesse beziehen. Die Erfahrung zeigt 
nun, daß insbesondere Hypothesen über den letzteren Sachver- 
halt sich immer wieder ’ in der Praxis als zu rigid erweisen, 
ganz abgesehen davon, daß unter den Konstrukteuren von Sai- 
sonbereinigungsverfahren kein Konsens über zu postulierende 
Traj ektor ieneigenschaf ten besteht und ein solcher auch nicht 
zu erwarten ist. Bei den in der statistischen Praxis domi- 
nierenden Verfahren wird nun versucht, diese Rigidität zu 
durchbrechen durch "Anleihen" von Konzepten aus Analysen, 
die im Frequenzbereich beheimatet sind. Es wurde oben schon 
dargelegt, warum eine solche Vorgehensweise grundsätzlich un- 
befriedigend ist. Will man derartigen "trial and error "- 
Prozeduren entgehen, dann scheint nur eine Analyse, die 
ausschließlich im Frequenzbereich durchgeführt wird, in Fra- 
ge zu kommen. Im Rahmen einer solchen Betrachtungsweise stellt 
sich dann die Saisonbereinigung als reines Filter-Problem dar. 
Dagegen lassen sich zwei Einwände erheben; 

Zunächst setzt sich eine derart vollzogene Saisonbereinigung 
möglicherweise der Kritik aus, ohne "ökonomisches Modell" zu 
operieren. Es ist jedoch leicht einzusehen, daß eine solche 
Kritik unberechtigt ist, da sie stillschweigend von der Vor- 
aussetzung ausgeht, bei den traditionellen Komponenten-Mo- 
dellen handele es sich tatsächlich um ökonomische Modelle. 
Genau besehen ist das jedoch keineswegs der Fall: den einzel- 
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nen Komponenten werden lediglich ökonomische Begriffe beige- 
legt/ aber es wird z. B. nicht aus der Wachstums- und Konjunk- 
turtheorie deduziert, daß die Trajektorien der glatten Kompo- 
nente etwa gerade ein Polynom dritter Ordnung ist. Genau das 
müßte aber der Fall sein, wenn der Anspruch gerechtfertigt 
sein soll, daß die traditionellen Komponentenzerlegung schon 
ein ökonomisches Modell darstelle. Daß mit den herkömmlichen 
Komponenten-Modellen ökonomisch nichts erklärt wird, folgt 
schon allein aus der Tatsache, daß die einzelnen Komponenten 
lediglich als Funktionen der Zeit angesehen werden, aber nicht 
als Funktionen ökonomischer Variablen. Zweifellos wäre der 
letztere Ansatz äußerst wünschenswert, jedoch ist nicht zu 
sehen, wie man aufgrund des heutigen Standes der einschlägi- 
gen ökonomischen Theorien mit solchen Kausalmodellen zu prak- 
tikablen und allgemein akzeptierten Verfahren kommen soll. 

Auch zukünftig ist für einen solchen Ansatz grundsätzlich 
Skepsis anzumelden, wenn man sich vergegenwärtigt, daß Zeit- 
reihenprobleme in der Physik, Elektrotechnik usw. heute grund- 
sätzlich mittels Filtertechniken gelöst werden, also weitge- 
hend auf "physikalische" Modelle verzichtet wird, obwohl in 
diesen Disziplinen weit entwickelte Theorien zur Verfügung 
stehen. Vergleicht man den dort erreichten Stand der Theori- 
enbildung mit dem in den Wirtschaftswissenschaften, dann er- 
scheint die geäußerte Skepsis leider mehr als gerechtfertigt. 

Der zweite Einwand ist technisch-statistischer Natur. Offen- 
sichtlich setzen Konzepte wie Transferfunktion, Spektrum, 
Relationen zwischen Filter-Input-Reihe und Filter-Output- 
Reihe usw. zunächst stationäre stochastische Prozesse voraus. 
Ökonomische Zeitreihen können aber nicht ohne weiteres als 
Realisationen derartiger Prozesse interpretiert werden. Das 
bedeutet aber noch lange nicht, daß diese Konzepte damit für 
die Analyse nicht-stationärer Prozesse gegenstandslos wären. 
Sie lassen sich nämlich sinnvoll auf derartige Prozesse 
übertragen, wenn durch geeignete Definitionen die Klasse der 
nicht-stationären Prozesse hinreichend eingeschränkt wird. 
Dadurch ergeben sich z. B. Konzepte wie "Durchschnittsspek- 
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tren" (Vgl. GRANGER, HATANAKA, 11) oder "evolutionäre Spektren" 
(PRIESTLEY, 29) für die sich die Wirkungen linearer Filter mit 
für praktische Zwecke hinreichender Genauigkeit wiederum an 
ihren Transferfunktionen bestimmen lassen. Am Rande sei er- 
wähnt, daß es übrigens im Rahmen der modernen Filtertheorie 
gleichgültig ist, ob eine Zeitreihe als Realisation eines 
stochastischen Prozesses oder als deterministische Funktion der 
Zeit interpretiert wird. Somit ist auch dieser zweite Einwand 
praktisch nicht von allzu großer Relevanz, insbesondere auch, 
wenn man bedenkt, daß eine sinnvolle harmonische Analyse 
nicht-stationärer Prozesse schon aus dem Grund möglich sein 
muß, weil realiter stationäre Prozesse nicht existieren und 
trotzdem sinnvolle Resultate mittels einer solchen Analyse 
in ganz verschiedenen Disziplinen wie der Biologie, der Geo- 
physik usw. erzielt wurden. (Vgl. dazu BRILLINGER, 5) 

Eine ausschließlich im Frequenzbereich durchgeführte Saison- 
bereinigung setzt natürlich auch eine f requentielle Defini- 
tion dessen voraus, was als saisonale Variation einer Zeit- 
reihe verstanden werden soll. Allgemein wird man zunächst 
als saisonale Variation einer Zeitreihe diejenigen Varia- 
tionen charakterisieren, die eine Konzentration spektraler 
Masse im Spektrum einer Zeitreihe "um" die saisonalen Fre- 
quenzen ^^,k=1,...,6 (bei Monatswerten) verursacht. Wie oben 
ausführlich dargestellt wurde, ist damit noch keine Defini- 
tion des interessierenden Phänomens gegeben, da nicht a priori 
feststeht, was unter "um" zu verstehen ist, d. h. also, die 
relevanten Bandbreiten (und damit Frequenzbänder) bleiben bei 
einer solchen allgemeinen Charakterisierung Undefiniert. Wie 
oben ebenfalls gezeigt wurde, erscheint eine objektive (d. h. 
also durch Verwendung statistischer Optimalitätskriterien) 
Bestimmung dieser Bandbreiten grundsätzlich als unmöglich, 
jedenfalls solange man sich Gütekriterien der statistischen 
Schätztheorie bedient. Da im Frequenzbereich "Saison" erst 
dann definiert ist, wenn diese Bandbreiten fest liegen, bleibt 
im Prinzip nichts anderes übrig, als diese vor der Filterkon- 
struktion subjektiv zu wählen. So wie bei einer Saisonberei- 




72 



nigung im Zeitbereich für die einzelnen Komponenten bestimmte 
Trajektorien ausgewählt werden, so müssen bei einer Saisonbe- 
reinigung im Frequenzbereich die relevanten Bandbreiten ge- 
wählt werden. Damit sind dann auch die jeweiligen oben auf ge- 
führten Transferfunktionen vollständig definiert. Allerdings 
ist ein solcher Ansatz wesentlich flexibler als eine a priori 
Definition von Komponenten-Traj ektorien , da mit vorgegebenen 
Bandbreiten eine Vielzahl von durchaus sehr unterschiedlich 
verlaufenden Komponenten-Prozessen vereinbar ist, denen nur 
gemeinsam ist, daß ihre jeweiligen Spektren in den relevanten 
durch die Wahl der Bandbreiten festgelegten, Frequenzbändern 
im wesentlichen konzentriert sind. Für eine sinnvolle subjek- 
tive Vergleichbarkeit verschiedener rein frequentiell konzi- 
pierter Saisonbereinigungsverfahren wäre natürlich eine Eini- 
gung über diese Bandbreiten äußerst wünschenswert, denn nur 
eine derartige "Normierung" würde gewährleisten, daß bei Ver- 
fahrensvergleichen nicht Größen verglichen werden, die a 
priori, d. h. also hier: ex definitione, unvergleichbar sind. 

Liegen die Frequenzbänder fest, dann ist eine Entscheidung 

darüber zu fällen, welche Transferfunktion realisiert werden 
1 ) 

soll. Das wird ganz entscheidend davon abhängen, welche 
Zielsetzung einer Saisonbereinigung zugrundeliegt. Zwei Ziel- 
setzungen erscheinen naheliegend: 

1. Elimination der saisonalen Variation einer Zeitreihe. 

2. Darstellung, d. h. Isolation, der saisonalen Variation 
einer Zeitreihe. 

Für die Zielsetzung 1 . ist eine Realisierung der obigen 
Transfer Funktionen II oder IV erforderlich (die Transfer- 
funktionen III und V sind hier nicht sinnvoll, da ja der Ana- 
lyse grundsätzlich kein Komponentenmodell zugrundeliegt) , 
während für die Zielsetzung 2. die Transferfunktion I in Fra- 
ge kommt. Das Resultat von Filter Operationen mit den korres- 
pondierenden Transferfunktionen II oder IV ist schlicht eine 



1) Es wird weiter unten noch ausgeführt, was im vorliegen- 
den Zusammenhang unter dieser "Realisierung" verstanden 
werden kann . 
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"saisonbereinigte" Reihe. Eine solchermaßen vollzogene Saison- 
bereinigung ist die denkbar "direkteste", da keine anderen 
als die relevanten Saisonbänder tangiert werden, keine iterati 
ven Prozeduren sowie Verknüpfungshypothesen erforderlich sind, 
wie bei den traditionellen Komponentenmodellen. Selbstver ständ 
lieh ist die bei diesen Komponentenmodellen übliche "subtile" 
Unterscheidung zwischen einer "saisonbereinigten" Zeitreihe 
und einer "glatten" Komponente dabei nicht möglich bzw. sinn- 
voll, da ein solches Komponentenmodell nicht vor liegt und da- 
mit auch keine "glatte Komponente" definiert ist. Am Rande sei 
nur vermerkt, daß sich gegen diese Unterscheidung, die mitt- 
lerweile zur Gewohnheit geworden ist, einiges einwenden läßt. 
Im Prnzip ist sie subjektiv, genauer gesagt: ihre verfahrens- 
abhängigen numerischen Differenzen sind letzten Endes defini- 
torisch bedingt. Doch kann darauf hier nicht weiter eingegan- 
gen werden. Es sei lediglich noch darauf hingewiesen, daß auch 
bei den traditionellen Komponentenmodellen die Zielsetzung 
"glatte Komponente" im Rahmen eines Saisonbereinigungsverfah- 
rens als nicht sehr naheliegend erscheint, da eine solche 
Komponente prinzipiell, d. h. unabhängig von einer Saisonbe- 
reinigung, mittels eines Tiefpaßfilters ermittelt werden kann, 
dessen Abschneidefrequenz unterhalb der ersten saisonalen 
Frequenz liegt. Somit ist nicht einzusehen, weshalb man diese 
Komponente eigentlich als Resultat eines Saisonbereinigungs- 
verfahrens erhalten zu müssen glaubt. 

Mit den beiden oben formulierten Zielsetzungen wird man aller- 
dings in der Praxis nicht ganz auskommen. Vielmehr sind dabei 
jeweils noch zwei Nebenbedingungen zu beachten. Die erste Ne- 
benbedingung ist die Forderung nach einer Saisonbereinigung 
bis zum "aktuellen Rand" einer Zeitreihe bzw. die Forderung 
nach der Isolation der saisonalen Variation einer Zeitreihe 
ebenfalls bis zum "aktuellen Rand". Diese Forderung impliziert 
daß die anzuwendenden Filter - im Idealfall - keine negativen 
Phasenverschiebungen induzieren. Die zweite Nebenbedingung 
ist die Forderung nach "Stabilität" der saisonbereinigten 
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Werte am aktuellen Rand bzw. die Forderung nach "Stabilität" 
der ausgefilterten saisonalen Variation am aktuellen Rand 
einer Zeitreihe. Diese Forderung bedeutet, daß sich die bis- 
herigen bereinigten bzw. ausgefilterten Werte beim Hinzukom-: 
men neuer Zeitreihenwerte - im Idealfall - nicht ändern. Über 
die unbestreitbare praktische Relevanz dieser beiden Ziel- 
setzungen (Vgl. dazu SCHÄIFFEI^ 33 ) braucht hier nicht dis- 
kutiert zu werden. Sie stehen jedoch in enger Beziehung mit 
dem zu wählenden Filtertyp. Dies soll nun etwas ausführlicher 
dargestellt werden. 

Grundsätzlich können die relevanten Transferfunktionen sowohl 
mittels rekursiver als auch mittels transversaler Filter rea- 
lisiert werden. 

Wie bereits dargelegt, haben beide Klassen von Filtern ihre 
spezifischen Vor- und Nachteile. So können z. B. mit rekur- 
siven Filtern niedriger Ordnung (also bei kleinem N und M) 
Transferfunktionen realisiert werden, die mit Transversalfil- 
tern nur mit sehr viel höherer Ordnung (also großem M) zu 
realisieren sind. Dagegen haben Transversalfilter grundsätz- 
lich lineare Phasenfunktionen (also auch Null-Phase) , wäh- 
rend rekursive Filter prinzipiell nicht-lineare Phasenfunk- 
tionen aufweisen (Vgl. KAISER, 16). 

Zur Berücksichtigung der oben skizzierten prospektiven Ziel- 
setzungen erscheinen zunächst Transversalfilter geeignet, da 
mit ihnen die einschlägigen Schätzungen ohne Phasenverschie- 
bung durchgeführt werden können. Allerdings ist das nur dann 
der Fall, wenn diese Filter in symmetrischer Form verwendet 
werden. Mit solchen Filtern ist jedoch z. B. eine Saisonberei- 
nigung bis zum aktuellen Rand nicht möglich. Obwohl rekursi- 
ve Filter nicht-lineare Phasenfunktionen aufweisen (also nie 
eine Null-Phase im gesamten Frequenzbereich [0,tt]) , sind 
darunter solche zu finden, deren Phasenfunktionen um den 
Wert Null schwanken, so daß die Phasenverschiebungen für prak- 
tisch alle Frequenzen kleiner sind als die der Analyse zugrun- 
deliegende Zeiteinheit (Monat, Quartal, Halbjahr) . Solche 
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Phasenverschiebungen wirken sich praktisch nicht aus und kön- 
nen deshalb vernachlässigt werden. Somit könnte mit solchen 
Filtern tatsächlich eine Analyse bis zum aktuellen Rand durch- 
geführt werden. Außerdem wäre die Stabilitätsforderung ideal 
erfüllt, da sich durch Hinzukommen neuer Zeitreihenwerte 
frühere Ergebnisse überhaupt nicht ändern. Für Zwecke der 
Saisonbereinigung (und auch isolierter Darstellung saisonaler 
Variation) sind deshalb prinzipiell rekursive Filter mit zu 
vernachlässigenden Phasenverschiebungen zu entwickeln. 

Es sei hier ausdrücklich vermerkt, daß der geschilderte Sach- 
verhalt kein spezifischer Nachteil von rein im Frequenzbe- 
reich konzipierten Verfahren ist. Auch bei den traditionallen 
Komponentenmodellen tritt diese Problematik in Erscheinung: 
z. B. werden sowohl beim Berliner als auch beim ASA-II-V er- 
fahren am Ende einer Reihe asymmetrische Stützbereiche für 
die gleitenden Regressionen verwendet, um mit der Analyse 
bis zum aktuellen Rand vorzustoßen. Tatsächlich implizieren 
diese asymmetrischen Stützbereiche jedoch negative Phasen- 
verschiebungen, so daß z. B. der letzte bereinigte Wert 
richtigerweise nicht dem letzten Zeitreihenwert, sondern dem 
Zeitreihenwert eines früheren Zeitpunktes zuzuordnen wäre. 
Somit ist das oben dargelegte Phasenproblem bei diesen Ver- 
fahren nur scheinbar gelöst und die Analyse wird in Wirk- 
lichkeit nicht bis zum letzten verfügbaren Zeitreihenwert 
durchgeführt . 

An dieser Stelle ist noch einmal auf die oben geführte Dis- 
kussion über die Wahl der einschlägigen Bandbreiten zurück- 
zukommen. Es wurde dargelegt, daß sich diese nicht "optimal" 
bestimmen lassen, sondern jeweils vom Verfahrenskonstrukteur 
festzulegen sind. Es ist nun durchaus denkbar, daß sich die 
a priori gewählten Bandbreiten nicht durchhalten lassen. Das 
kann verschiedene Gründe haben. Beispielsweise können diese 
sich ungünstig auf den verwendeten Approximationsalgorithmus 
auswirken, also zu numerischen Instabilitäten führen, oder 
etwa zu nicht zu akzeptierenden Phasenfunktionen. Gewählte 
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Bandbreiten und Filterkonstruktionen können also durchaus kol- 
lidieren. In derartigen Fällen ist dann eine Revision der 
ursprünglich gewählten Bandbreiten unumgänglich, eventuell so 
lange, bis diese unerwünschten Nebeneffekte völlig verschwun- 
den sind. Die aus einem solchen Revisionsprozeß schließlich 
resultierenden Bandbreiten können dann als "optimal" in die- 
sem Sinne bezeichnet werden. Allerdings handelt es sich dabei 
nicht mehr um ein übliches Optimalitätskriterium im Sinne der 
statistischen Schätztheorie , sondern mehr um eine pragmatisch 
orientierte Optimalität. Für die oben diskutierte grund- 
sätzliche Definitionsproblematik der mittels Filter numerisch 
zu konkretisierenden Phänomene würde das implizieren, daß die- 
se Phänomene im Frequenzbereich prinzipiell so definiert wer- 
den, daß diese Definitionen filtertheoretisch operational sind, 
also zu praktikablen Filtern führen. 

Die vorstehenden Untersuchungen zeigen, daß auch eine "Saison- 
bereinigung" mit den Werkzeugen der modernen Filtertheorie 
recht subtile Überlegungen erfordert. Insbesondere dürfte 
klar geworden sein, daß auch bei diesem Zugang "Freiheitsgra- 
de" für den Verfahrenskonstrukteur vorhanden sind, als deren 
wichtigster wohl die Wahl des Approximationskriteriums anzu- 
sehen ist. Jedoch ist eine derartige Analyse auf allen Stu- 
fen klar durchschaubar und theoretisch fundiert. Das hat zur 
Konsequenz, daß numerische Differenzen zwischen verschiedenen 
rein filtertheoretisch konzipierten Verfahren relativ leicht 
zu "orten" sind: verschiedene Wahl der relevanten Frequenz- 
bänder, verschiedene zugrundegelegte Transferfunktionen (und 
damit auch unterschiedliche Zielsetzungen) , verschiedene Fil- 
tertypen, verschiedene Approximationskriterien sowie verschie- 
dene Implementierungen kommen dafür alternativ oder in Kombi- 
nation in Frage. Außerdem kommen solche Verfahren mit einem 
Minimum an Hypothesen aus. Es ist zu erwarten, daß der filter- 
theoretisch orientierte Zugang zum Problemkreis "Saisonberei- 
nigung" zu einer Versachlichung und vor allem zu einer Ver- 
einheitlichung der heute zuweilen recht konfus anmutenden 
Diskussionen führen dürfte. 
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Die einzelnen Phasen einer "Saisonbereinigung" im diskutierten 
Sinn, die dem Verfasser als zu erstrebendes Ideal vor schwebt, 
seien der Übersichtlichkeit halber noch einmal in folgendem 
Ablaufdiagramm schematisch dargestellt: 
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3.4. Ein neues Verfahren zur 

Saisonbereinigung ökonomischer 
Zeitreihen 



Entsprechend den obigen Ausführungen wird nun hier ein Ver- 
fahren skizziert, dessen Ziel es ist, eine Filterfolge zu 
finden, welche die vorher erwähnte "ideale" Transf erfunktion 
IV realisiert. Die Zielsetzung eines solchen Verfahrens kann 
charakterisiert werden als "Elimination der saisonalen Va- 
riation einer Zeitreihe". Aus praktischen Erwägungen heraus 
sollten dabei jedoch zwei Nebenbedingungen beachtet werden, 
wie schon oben angeführt, nämlich: 

1) Die Analyse ist möglichst bis zum aktuellen Rand einer 
Zeitreihe vorzunehmen, d. h. Phasenverschiebungen sind 
tunlichst zu vermeiden. Diese Nebenbedingung soll kurz, 
wenn auch nicht ganz zutreffend, als "Nullphase" be- 
zeichnet werden. 

2) Bereits gefilterte Werte sollen sich beim Hinzukommen 
neuer Originalwerte nicht verändern. Diese Nebenbedin- 
gung sei kurz "Randstabilität" genannt. 

Die letztere Bedingung erfordert zwingend den Einsatz 
von rekursiven Filtern. 



Die gesuchte Transferfunktion ergibt sich nun durch eine 
sechsfache Serienschaltung des Kerbenfilters (für Monats- 
daten) . Man erhält dann: 



Ts(z) 



irb, ( z- e^^ k) ^ ( z- e ^^k)^ 
k=1 ^ 

6 3 



TT TT 

k=1 i=1 



(z-Pi,^) (z-p. 



iX 

/ z = e 



|A |< TT 



mit 



bk : = Normierungskonstante des k-ten 
Saisonfilters 

= i-ter Pol des k-ten Saisonfilters 
i= 1,2,3 k=1,2,...6 
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Der gesamte Filterprozeß kann schematisch v/ie folgt darge- 
stellt werden: 



X (t)-. 


T 

"1 


— »y-, (t)— . 


^2 


— ^^) 







♦VgCt) 



Die Zeitreihe x(t) wird also sukzessive durch sechs Kerben- 
filter geschickt mit den Transferfunktionen T , . . . ,T , deren 

27Tk ^ 

Nullstellen bei A = ^ ~ ^ liegen. Dabei ist 

jede Transferfunktion T (k = 1, 2, ... 6) von der Ordnung 

(M = 6, N = 6), saisonbereinigte Zeitreihe. 

Die i-te Filtergleichung des k-ten Saisonfilters kann dabei 
durch folgende Überlegung abgeleitet werden: 

Es ist 

Tg (z) = - e"\)(z-e-""k). 

(- - Pik) (^ - ?ik) 

-2 - 1 

z -2z cos A, + 1 
_ k 

IPikl^ - (p.j^ + p.^) z'”' + 1 



Da die Beziehung gilt: 



V 

‘ik 



(z) 



mit 




(z) (z) 



i — 1/2,3, k — 1,2, ...6 



Y. (z) : = z -Transformierte des Filteroutputs der 

1 K 

i-ten Filtergleichung des k-ten 
Saif onf ilter s 



X., (z) : = z-Transf ormierte des Filterinputs der 

1 K. 

i-ten Filtergleichung des k-ten 
Saisonfilters 
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folgt: 



Yik(2)UPikl^ - (Pik + Pik» 



= X^j^(z) [z ^ - 2 z ^ cosAk + 1] 



und nach Rücktransformation in den Zeitbereich: 



Yik(t) = x^j^(t) - 2 cosAj^ ^ik^^"^^ ■*" ^ik^^“^^ 

+ (Pik + Pik> ^ik^^"^^ - iPiki^ 
i = 1 ,2,3 , k = 1 ,2, . . .6 

(Vgl. zur letzteren Formel CADZOW, 6 ) 

Die Abb. 34, 35 zeigen die Transferfunktionen |t^(A) | für 
ot = 72^ und A = 0,01 , A = 0,1 . 

Zweifellos entspricht diese Transferfunktion noch nicht den 
theoretischen Vorstellungen, da auch sie, wie schon der ein- 
fache Kerbenfilter, "Nasen" aufweist. Die Schritte zur Besei 
tigung dieser Erscheinung sollen hier nicht dargestellt v/er- 
den, genau so wenig wie die recht umfangreichen und teilv/ei- 
se komplizierten Untersuchungen zu Fragen der numerischen 
Stabilität sowie der Phaseneigenschaften und den "Einschwing 
Vorgängen" des Verfahrens. Alle diese Probleme werden aus- 
führlich in einer gesonderten Publikation des Verfassers be- 
handelt. Darin finden sich auch Anwendungsbeispiele. 
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3.5. Ein neues Verfahren zur 

Schätzung einer "glatten Reihe" 

Die beim vorherigen Verfahren intendierte Zielsetzung einer 
"Elimination der saisonalen Variation einer Zeitreihe" ist 
nicht die einzig denkbare bzv/. sinnvolle Zielsetzung für 
derartige Bereinigungsverfahren. Häufig interessiert sich 
vor allem der Praktiker in erster Linie für diejenigen In- 
formationen einer Zeitreihe, die dem Niederfrequenzbereich 
zuzuordnen sind. Eine Ausfilterung solcher Informationen 
erfordert natürlich eine völlig andere Filterkonstellation 
als das vorher skizzierte Saisonbereinigungsverfahren. Wie 
unter 3.3 begründet vnarde, ist dazu die Konstruktion eines 
Tiefpaßfilters notwendig, dessen Abschneidefrequenz unter- 
halb der ersten Saisonfrequenz liegt. 

Man könnte nun versuchen, einen solchen Filter als optima- 
len IIR-Filter zu konstruieren. Dabei ergeben sich aller- 
dings zwei Probleme. Bei niedriger Filterordnung macht sich 
je nach vorliegender Zeitreihe die Differenz |t( 0)-1 | un- 
angenehm bemerkbar, d. h. die geschätzte "glatte Reihe" 
weist einen bias auf. Wird deshalb nun die Filterordnung 
erhöht, dann reduziert sich der bias, aber gleichzeitig 
verschlechtern sich die Phaseneigenschaften, so daß man 
keine phasenfreien Schätzungen erhält (vgl. Abb. 36, 37). 
Die Nebenbedingung "Nullphase" sollte aber hier aus den 
selben Gründen wie beim obigen Saisonbereinigungsverfahren 
erfüllt sein. 

Die Phaseneigenschaften lassen sich grundsätzlich verbes- 
sern, wenn der Übergang zwischen Durchlass- und Stoppband 
nicht abrupt erfolgt, sonder "glatter", was z. B. dann der 
Fall ist, vjenn ein Übergangsband definiert wird, in dem die 
Transferfunktion monoton (z. B. linear) abnimmt. Wie ein- 
schlägige Untersuchungen zeigten, ist die Phasenverschie- 
bung zv^ischen Originalreihe und ausgefilterten Reihen 
zwar geringer als vorher, jedoch ist diese Verbesserung 
nur eine unwesentliche. 




82 



Weitere Versuche in dieser Richtung brachten leider keine 
durchweg befriedigenden Resultate. 

Eine andere Möglichkeit, zu einem brauchbaren Tiefpaßfilter 
zu kommen, besteht darin, einen aus der Klasse der "maximal 
flachen" Butterworth-Tief paßfilter auszuwählen. Diese Fil- 
ter haben die Eigenschaft, daß alle Ableitungen der Trans- 
ferfunktion bei X = 0 den Wert Null haben (vgl. RABINER - 
GOLD, 30) . Außerdem verfügen sie über eine Normierungskon- 
stante, die immer so gewählt werden kann, daß T (0) =1 ist. 
Somit tritt bei Verwendung eines solchen Filters das vorher 
erwähnte bias-Problem nicht auf. Die Transferfunktion eines 
solchen Filters lautet z. B.: 



k 

TT ( z - p ) 

i=1 ^ 



mit 



b 



W (1 - p.) 

i=1 ^ 



SO daß gilt T(1) = 1, d. h. T(A) ist Eins bei X = 0. 

(Vgl. CADZOW, 6 ) . Da aber auch bei diesem Filter der Durch- 
lassbereich relativ klein ist (vgl. Abb. 38), ergeben sich 
unerwünschte Phasenverschiebungen, v;ie die Abb. 39 zeigt. 

Es bleibt also das Problem zu lösen, einen Tiefpaßfilter zu 
konstruieren, der sowohl Nullphase hat als auch die Bedin- 
gung T(X) = 1 für X = 0 erfüllt und außerdem im Durchlassbe- 
reich völlig flach verläuft. Die Lösung dieses Problems liegt 
in einer geschickten Kombination von unter Umständen mehre- 
ren rekursiven Filtern. Sie soll hier ebenfalls nicht darge- 
stellt werden, da hierüber eine spezielle Arbeit des Verfas- 
sers in Vorbereitung ist. 
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3.6. Rekursive Filter und 
Zeitreihenprognose 

Oft wird statt dem Ausdruck "Prognoseverfahren" auch der Aus- 
druck "Prognosefilter" verwendet Um etv/aige Verv/echselun- 

gen zu vermeiden, sei gleich darauf hingewiesen, daß die oben 
behandelten Filter natürlich reine Analyseinstrumente darstel- 
len und inhaltlich mit irgendv/elchen Prognoseverfahren über- 
haupt nichts zu tun haben. Wenn in diesem Abschnitt rekursive 
Filter mit Prognose in Verbindung gebracht werden, dann nicht 
deswegen, weil diese möglicherweise an sich schon für prog- 
nostische Zwecke umfunktioniert werden könnten, sondern aus 
dem Grund, vzeil es m. E. möglich und sinnvoll sein dürfte, 
diese mit existierenden Prognose-Verfahren zu verknüpfen. 

Speziell v/ird dabei hier an Prognosemodelle vom Typ BOX- 
JENKINS gedacht. Bei den univariaten Modellen V7ird dabei 
postuliert, daß sich eine Zeitreihe X^, die im Allgemeinen 
eine Realisation eines nicht-stationären stochastischen Pro- 
zesses ist, wie folgt darstellen läßt: 





(p(B) 


= e(B) 


mit 


g)(B) 


: = 0(B) (l-B)*^ 




0(B) 


: = 1 - 0^B - 0^B^ 0pB> 




e(B) 


: = 1 - Ö-B - e„B^ - - - ö b‘ 

1 Z q 



wobei B ein Verschiebungsoperator ist, für den gilt: 

X. 






t-j 



a^ ist ein reiner Zufallsprozeß ("weißes Rauschen"). Diese 
Darstellung impliziert die Vorstellung, daß eine Zeitreihe 
nac h damaliger Differenzenbildung als autoregressiver-glei- 
tender Durchschnittsprozeß darstellbar ist (ARIMA-Prozeß 

1) Dies hängt wohl damit zusammen, daß die Prognoseverfah- 
ren von WIENER für stationäre Prozesse auch unter dem 
Namen WIENER-Filter bekannt sind. 




84 



der Ordnung (p, d, q) ) . Dieser univariate Ansatz läßt sich 
leicht erweitern. Liegen z. B. zwei Reihen Y^) vor, dann 

läßt sich möglicherweise dieser bivariate Prozeß durch ein 
"Transf erfunktionen-Modell" der Ordnung (r, s) dar stellen: 

6(B) = o)(B) 

mit 

6 (B) : = 1 - 6^B 6^b'^ 

g 

0)(B) : = o)q - (jü^B - - - oü^B 

wobei b eine Phasenverschiebung zwischen beiden Prozessen 
bedeutet . 

Praktisch wird obige Relation zwischen Y^ und nicht un- 
gestört gelten, sondern von einem "Rausch"-Prozeß über- 
lagert sein, der selbst als ARIMA-Prozeß darstellbar ist, 
so daß also gilt: 

6 (B) Y^ = cü(B)X^_j^ 4- 

(B) 0(B) 

(Vgl. BOX-JENKINS, 4) 

Ein interessanter Aspekt dieser Transferfunktionenmodelle 
liegt nun darin, daß sie es gestatten, zur Prognose einer 
Reihe Informationen über lead-lag-Beziehungen mit anderen 
Reihen heranzuziehen, was für die Prognose ökonomischer Zeit- 
reihen von besonderer Bedeutung sein dürfte. 

Viel diskutiert wurde in den letzten Jahren das Problem der 
Konstruktion konjunktureller Indikatoren, sowohl zu Diagno- 
se- als auch zu Prognosezwecken. Wie oben angedeutet vTurde, 
ist es möglich, mit Hilfe einer speziellen Kombination von 
rekursiven Filtern die konjunkturelle Komponente von Zeit- 




85 



reihen isoliert darzustellen. Eventuell vorhandene Phasen- 
verschiebungen zwischen diesen Komponenten dürften sich mit 
Hilfe der Kreuz Spektralanalyse ermitteln lassen. Es ist zu 
erwarten, daß die sonst beim Einsatz der Spektralanalyse für 
ökonomische Reihen festzustellenden Schwierigkeiten dabei 
nicht auftreten, z. B. gibt es kein " leakage" -Problem im 
Niederfrequenzbereich das von einer Trendkomponente her- 
rührt. Die Phasenbeziehungen zwischen einzelnen Konjunktur- 
Komponenten dürften sich daher ziemlich problemlos schät- 
zen lassen. Kann man so "führende" Konjunkturreihen er- 
mitteln, dann müßte es auch möglich sein, damit sinnvolle 
Transferfunktionenmodelle unter Beachtung eben dieser Pha- 
seninformationen für die isolierte konjunkturelle Variation 
wichtiger Ökonomischer Zeitreihen zu formulieren und da- 
mit die Prognostizierbarkeit konjunktureller Prozesse mög- 
licherweise entscheidend zu verbessern. 

Es soll hier nicht verschwiegen werden, daß dieser Problem- 
kreis aus Zeitgründen im Rahmen dieses Forschungsprojektes 
nicht oder nur fragmentarisch bearbeitet werden konnte. 

Hier liegt ein weites Feld zukünftiger Aktivitäten, auf 
dem sich Filtertheorie und Prognosetheorie verbinden. 
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Abb , 1 ; Rechteckfilter, N=23, fj^=0,05 (mit Unterschätzung des Trends) 



<196. 




ftbb. 2: Rechteckfilter, N=23, f;^=0,05 (ohne Unterschätzung des Trends) 
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Abb. 4: Kaiser-Filter 
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Abb. 5; Geschätzte glatte Komponente des Indexes der industriellen Nettoproduktion, gesamte Industrie 
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G 0 schätzte glatte Komponente des Indexes der industriellen Nettoproduktion ^ gesamte Industrie 
(Dez. 1962 - Jan. 1975), Kaiser-Filter, N=23, f_=0,05 
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Gleitender Durchschnitts-Filter 
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Abb . 15: Transferfunktion des optimalen FIR-Filters der Fil- 
terlänge 99 und der Gewichtsrelation 10 : 1 
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Nettoprodviktion verarbeitende Industrie 
gefiltert mit FIR-Tiefpaß (Filterlänge 99) 
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Abb . 21; o Simulationsreihe, A gefiltert mit IIR-Tiefpaß, 
f =0.05, n=l 
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22; o Simulationsreihe, A gefiltert mit IIR-Tiefpaß, 
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23: o Simulationsreihe, A gefiltert mit IIR-Tiefpaß, 
f =0.05, n=45 
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26; o Simulationsreihe, A gefiltert mit IIR-Bandfilter 
n=l, [0.02, 0.04] 
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Abb. 27; o Simulationsreihe, A gefiltert mit IIR-Bandfilter 
n=4,' [0.02, 0.04] 
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28: o Simulationsreihe, A gefiltert mit IIR-Bandfilter 
n=l , [0.025, 0.035] 
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:ionsreihe, A gefilteirt mit IIR-Bandf ilter 
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30; o Simulationsreihe, A gefiltert mit IIR-Bandfilter 
n=l, [0.029, 0.031] 



3094 



-115- 




31 : o Simulationsreihe, A gefiltert mit IIR-Bandf ilter 
n=4, [0.029, 0.031] 
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37; o Arbeitslose in der Bundesrepublik 
A gefiltert mit obigem Tiefpaßfilter 
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38; Transferfunktion Butterworth-Tiefpaßfilter 
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39 : o Arbeitslose in der Bundesrepublik 
A gefiltert mit obigem Tiefpaßfilter 
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